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Abstract: If dark matter couples directly to a light force mediator, then it may form bound 
states in the early universe and in the non-relativistic environment of haloes today. In this work, 
we establish a field-theoretic framework for the computation of bound-state formation cross-sections, 
de-excitation and decay rates, in theories with long-range interactions. Using this formalism, we 
carry out specific computations for scalar particles interacting either via a light scalar or vector 
mediator. At low relative velocities of the interacting particles, the formation of bound states 
is enhanced by the Sommerfeld effect. For particle-antiparticle pairs, we show that bound-state 
formation can be faster than annihilation into radiation in the regime where the Sommerfeld effect 
is important. The field-theoretic formalism outlined here can be generalised to compute bound- 
state formation cross-sections in a variety of theories, including theories featuring non-Abelian 
(albeit non-confining) interactions, such as the electroweak interactions. 
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1 Introduction 


Dark matter (DM) with long-range interactions, mediated by a light or massless force carrier, 
appears in diverse theories motivated on various grounds. Let us mention a few important examples. 
Dark matter with sizable self-interactions, mediated by a light particle, can explain the observed 
galactic structure better than collisionless DM [1-10]. Asymmetric DM [11-14] - motivated in part 
by the similarity of the dark and the ordinary matter abundances - resides, in most implementations, 
in a hidden sector that often includes long-range interactions [15-20]. Dark matter which dissipates 
energy via its coupling to a light force mediator may provide a dynamical explanation for some of 
the observed scaling relations governing haloes [21-24], and other features [25, 26]. Inside haloes, 
the inelastic scatterings of either symmetric or asymmetric DM with long-range interactions can 
produce radiative signals [27-35], which can potentially account for anomalous excesses observed in 
the radiation backgrounds [27-32]. Moreover, long-range DM-nucleon scattering implies a different 
interpretation of the direct-detection data than the commonly assumed short-range scattering [36- 
39]. Notably, the long-range character of DM interactions is relevant not only for theories involving 
hidden sectors; even the electroweak interactions of the Standard Model manifest as long-range if 
DM is heavier than a few TeV [40-42]. Clearly, long-range interactions play a central role in the 
venture to identify DM. In order to extract accurate predictions for the DM phenomenology, it is 
then essential to fully understand their implications. 

Long-range interactions typically imply the existence of bound states. The formation of DM 
bound states in the early universe and/or in the dense environment of haloes today affects the 
phenomenology of DM in many important ways. In the early universe, symmetric DM may form 
unstable bound states, whose decay contributes to the DM annihilation rate and affects the DM 
relic abundance [43]. Asymmetric DM may form stable bound states [15-20, 34, 44-52], which affect 
all manifestations of DM today, including the DM self-scattering in haloes [19, 46-48], the expected 
indirect-detection [29-34] and direct-detection [53] signals, as well as the kinetic decoupling of DM 
from dark radiation [54]. Inside haloes, DM bound states - whether they are stable or unstable 
- may form radiatively and yield detectable signatures [28-30, 35]. Radiative signals may also be 
produced in transitions between the bound-state energy levels [31-33], or other related inelastic 
processes [34]. Moreover, the formation of DM bound states may be detectable at the LHC [55]. 
These rich phenomenological implications strongly suggest that it is critical to accurately account 
for the formation of DM bound states, in order to constrain the DM properties and eventually 
detect DM. This work is a step toward this goal. 

There are, of course, two different classes of bound states: Those arising due to non-confining 
interactions, such as the atomic bound states in QED, and those arising due to confining inter¬ 
actions, such as the hadrons in QCD. Particles charged under a conhning force always combine 
into hadronic states, roughly once the kinetic energy in their center-of-mass frame (or the tem¬ 
perature of their plasma) drops below the confinement scale. On the other hand, in the case of 
non-confining interactions, the efficiency of bound-state formation (BSF) depends on the corre¬ 
sponding cross-sections and the details of the thermodynamic environment. Here we shall consider 
bound states due to non-confining interactions, and calculate the cross-sections for their formation 
in the non-relativistic regime, which is relevant for cosmology and DM indirect detection signals.^ 

The formation of DM bound states in the early universe and inside the non-relativistic envi¬ 
ronment of haloes differs from BSF in colliders in some important ways. In high-energy colliders, 
the initial-state particles are highly relativistic.^ However, BSF is more efficient when the relative 
velocity of the interacting particles is lower than the expectation valne of the relative velocity of the 
particles inside the bound state. Equivalently, this is when the kinetic energy in the center-of-mass 

^Another class of bound states - non-topological solitons - has also been considered in the context of DM [2, 56—60]. 
^See, however, Ref. [61]. 
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(CM) frame is lower than the binding energy; clearly, this lies within the non-relativistic regime. 
In this regime, the long-range interaction distorts the wavefunctions of the incoming particles, 
which cannot be approximated by momentum eigenstates (plane waves). This is the well-known 
Sommerfeld effect [62]. If the interaction is attractive, the Sommerfeld effect enhances the cross- 
section for any process the two particles may participate in, including the formation of bound states. 
It follows that, cosmologically, DM bound states form most efficiently after the temperature of the 
dark plasma drops below the binding energy. While the formation of bound states in the early 
universe eventually freezes out due to the cosmological expansion, bound states may again form 
efficiently in today’s dense and non-relativistic haloes. In either case, to accurately estimate the for¬ 
mation of DM bound states, we must account for the Sommerfeld effect, which is a non-perturbative 
phenomenon, as is of course the very existence of bound states. 

In this paper, we establish a field-theoretic framework for the calculation of BSF cross-sections 
and other level-transition rates, as well as the decay rates of unstable bound states. Then, we carry 
out computations for specific interactions. We organise our work as follows. 

We begin, in Sec. 2, with reviewing how to determine the wavefunctions of the two-particle 
states and the bound states in the presence of a long-range interaction. We derive the Bethe-Salpeter 
equation for the wavefunctions, and reduce it to the Schrodinger equation using the instantaneous 
approximation in the non-relativistic regime. In Sec. 3, we determine the amplitudes for transitions 
between energy levels with the emission of a force mediator; this includes the radiative capture 
to a bound state. In the fully relativistic regime, we express the amplitudes for such processes in 
terms of the Bethe-Salpeter wavefunctions describing the initial and final states, and a perturbative 
interaction. We then employ the instantaneous and the non-relativistic approximations, to express 
the transition amplitudes in terms of the Schrodinger wavefunctions. In Sec. 4, we repeat the 
analysis for the decay of unstable bound states into radiation, as well as for the (co-)annihilation 
of unbound pairs of particles into radiation - two closely related processes. While we lay out our 
formalism in terms of scalar particles, it is straightforward to extend it to include fermionic species. 

We continue by applying our formalism to specific interactions. In Sec. 5, we consider scalar 
particles interacting either via a scalar or a vector mediator, and calculate the cross-sections for 
the dominant radiative capture to a bound state. We estimate the range of validity of our com¬ 
putations, using the unitarity bound on the inelastic cross-section. For bound states made of 
particle-antiparticle pairs or pairs of self-conjugate particles, we compare BSF with annihilation, 
and show that BSF can be the dominant inelastic process in the regime where the Sommerfeld 
effect is important; we sketch this comparison in Fig. 1. In addition, we calculate the decay rates of 
particle-antiparticle bound states into force mediators. We cast our results in terms of a minimal 
parametrisation, which makes their potential implications more transparent, and summarise them 
in table 1. We conclude with a discussion of the phenomenological implications of DM bound-state 
formation in Sec. 6, and present many of the detailed calculations in the appendices. 

The field-theoretic formalism developed in this work has several advantages in comparison to 
previous quantum mechanical calculations [29, 35, 51, 63, 64]. It can accommodate the possibility 
of DM coupled to non-Abelian interactions. Such interactions can convert the incoming parti¬ 
cles into different species which may subsequently form bound states; importantly, DM coupled 
to the electroweak interactions of the Standard Model belongs to this category. Moreover, the 
field-theoretic approach allows for a systematic inclusion of higher-order corrections, both in the 
interaction strength and in the momentum transfer between the interacting degrees of freedom. 

We emphasise that BSF is distinct from processes such as the direct annihilation into mediators 
or elastic scattering, in which particles coupled to a long-range interaction may participate. While 
all these processes are influenced by the Sommerfeld effect, the final-state particles are obviously 
different. Field-theoretic treatments of the annihilation processes, analogous to the formalism 
presented here for BSF and discrete level transitions, have been presented in Refs. [65, 66]. 
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Table 1. Summary of the bound-state-formation and annihilation cross-sections, the decay and de¬ 
excitation rates, computed in Sec. 5 , for scalar particles interacting via a light scalar or vector media¬ 
tor. The annihilation cross-sections and the rates of decay into (^-mediators refer to unbound and bound 
particle-antiparticle pairs respectively; all other formulae apply to any pair of particles. The cross-sections 
are normalised to CTc = {-kcP/ fj?) x 27 r^/(l — where a is the hue-structure constant entering the 

Coulomb potential, /r = mi m2/(mi -f m2) is the reduced mass of the interacting species, and = a/urei, 
with Urei being the relative velocity of the incoming particles. Also, ?7i_2 = mi,2/(mi -|-m2), and <71,2, Ci,2 are 
the couplings of the interacting particles to the force mediators, as described by the Lagrangian densities 
of Sec. 5 . In our computations, we have neglected the mediator mass. Unitarity suggests that the range of 
validity of the above computations is a < 0 . 5 . 


2 Bound-state and two-particle state wavefnnctions 

We shall consider two scalar particle species xi and X 2 , interacting via a light or massless force 
mediator (p. Specific interaction Lagrangians will be introduced in Sec. 5 (c.f. Eqs. (5.13), (5.14) 
and (5.48)). In this section, we aim to determine the wavefunctions which describe two-particle 
states and bound states with the quantum numbers of xi and X 2 - Our presentation draws largely 
from the pedagogical discussions of Refs. [67, 68]. 

Let us first introduce some notation. We denote the masses of xi, X 2 and tp by m-i, m 2 and 
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Figure 1. Comparison of the cross-sections for annihilation (dashed blue lines) and the dominant radiative 
capture process to a bound state (solid purple lines), for a scalar particle-antiparticle pair. The cross- 
sections times relative velocity, aVrei, are normalised to the perturbative value for s-wave annihilation, cto- 
Left: For a scalar mediator, ao = 'ko? f[P'. The leading capture process is to the excited state {210} and 
consists dominantly of the J = 2 partial wave. At ^ 3> 1, the ratio of the bound-state formation and 
annihilation cross-sections is o-ggp/cTann — 0.066. Right: For a vector mediator, ao = 'kcP The 

leading capture process is to the ground state {100} and consists dominantly of the J = 0 and J = 2 partial 
waves. At ^ 1, O-ggp/Uann — 1.56. 


m^p respectively. We define the total and the reduced xi, X’l masses, 

m = mi -|- m2 , 

TOi m2 

R = -^- 

TOi -I- m2 

Obviously, /i ^ m/4. In the following, stands for a xi — X 2 bound state, of total momentum 

Q and energy -|- M^, where n denotes collectively all the discrete quantum numbers 

characterizing the bound state, and Mn < m is the bound-state mass. I^q.q) stands for a xi ~ X 2 
unbound two-particle state, with total momentum Q, expectation value of relative velocity v^ei = 
q/ fi, and energy ^ m. (As is common in scattering theory, we shall often refer to the unbound 
states as scattering states.) Moreover, |i^q) stands for a ip particle state with momentum Q and 
energy +m'^. 

2.1 The Bethe-Salpeter wavefunctions 

We shall now introduce the Bethe-Salpeter wavefunctions that will appear in the fully relativis¬ 
tic version of the transition amplitudes of sections 3 and 4. The Bethe-Salpeter wavefunctions 
are related to the more familiar Schrddinger wavefunctions, which we will use in evaluating the 
amplitudes of interest in the non-relativistic regime. 

We are interested in the following wavefunctions: 


( 2 . 1 ) 

( 2 . 2 ) 


and 


^'q,n(a;i,a:2) = (n|Txi(xi)x2(a;2)|Sq.n) , (2.3) 

^q,n(a;i,a;2) = {BQ,n\Txli.Xi)xUx2)\^) (2.4) 

®q.q(a;i,a;2) = {^\Txi{xi)x2ix2)\^Q,q) , (2.5) 

®q,q(a^Ta;2) = {UQ.q\Tx\ixi)xUx2)\^) , (2.6) 
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where T is the time-ordering operator, and |ri) is the vacuum of the interacting theory. If Xi,X 2 
are a particle-antiparticle pair, in the above definitions we replace xi Xi X 2 —t X^- Note that * 
does not denote complex conjugation, for which we shall use the symbol *, as usual. In fact 

= {n\fxi{Xl)x2{x2)\Bct,nY , (2.7) 

= (f 2 |fxi(a^i)x 2 (x 2 )|Z^Q,q)* , ( 2 . 8 ) 

where T is the anti-time-ordering operator. 

We define the coordinate transformation and its inverse 


X = xi - X2, X = 771 X 1 -I- 772 X 2 , (2.9) 

xi=X-|- 772X, X2 = X — rjix, (2.10) 


where 771 -|- 772 = 1 for the Jacobian to be I, and we choose specifically 

mi^2 

^— . 
nil + 7772 


( 2 . 11 ) 


In the non-relativistic regime, this choice will enable us to separate the motion of the CM from the 
relative motion. Using the 4-momentum operator P, we obtain 


Xi(xi) = exp(jPX)xi( 772 x) exp(- 7 PX) , (2.12) 

X2 {x 2) = exp{iPx)x2{-Vix) exp{-iPX) . (2.13) 


Then, the wavefunction of Eqs. (2.3) becomes 

^'Q,n(a;i,X 2 ) = 6»(x°)(ll|xi(a7i)x2(a;2)|^Q.n) + 0{-x°){n\x2ix2)xiixi)\BQ,n) 

= 6'(x°)(fl|e*^^xi(»72a;)x2(-?7ia;)e"*^^|PQ,„) -f 6'(-a;°)(ll|e*^^X2(-?7ia:)xi(»72a;)e"*^^|PQ,„) 

= e"*‘^^(n|rxi(772x)x2(-?7ia^)I^Q.n) 

= e-*<3^TQ,„(x) , (2.14) 

where it is understood that = tuq^n, and we defined 

^'Q,n(a;) = (n|rxi(772x)x2(-??ia;)|PQ,n) . (2.15) 


This is the first step in separating the motion of the CM from the relative motion. For notational 
simplicity, we are using the same symbol for 4 >q_„(xi,X 2 ) and 'I'q^„(x). We also define the Fourier 
transforms 

^Q.n(a:) = J ^ ^qAp) , ^Q,n(p) = J <Px vI/q.„(x) P?'" . ( 2.16) 

We repeat the above for the amplitudes of Fqs. (2.4) - (2.6), and summarise the definitions in 
appendix A. 


2.2 The 4-point Green’s function and Dyson-Schwinger equation 

Consider the 4-point Green’s function 

GN)(xi,X2,7/i, 7/2) = (ll|Txi(xi)x2(a;2)Xi(j/i)xk2/2)|ll) , (2.17) 

and let lU(xi,X 2 , 771 , 7 / 2 ) be the perturbative 4-point interaction kernel between xi and X 2 . Then, 
GN) satisfies the Dyson-Schwinger equation, 

G(^^(xi,X2,7/i,7/2) = -5'i(xi - yi)S 2 {x 2 - 7 / 2 ) + 

J (Pzi (Pz'i Ipz2 (Pz2 5'i(xi - Zi)S' 2 (x 2 - Z2)W{zi, Z2; z'l, z'^) Si{z[ - 7/l)S'2(4 “ ^ 2 ) + . . . , 

(2.18) 
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Figure 2. Diagrammatic representation of the Dyson-Schwinger equation (2.19) for the 4-point function 
G^^\xi,X2,yi,y2)- —9 — and —•— stand for the Xi and X 2 full propagators, respectively. 


where 81,82 are the full propagators for Xi)X 2 - Symbolically, the above series can be written as 

G(4) = S 182 + 8182 W 8182 + 8182 W 8182 W 8182 + ■■■ 

= 8182 + 8182 W {8182 + 8182 W 8182 + . . .) 

= 8182 + 8 i 82 WG^*'> . (2.19) 

Equation (2.19) is sketched in Fig. 2. 

Due to translational invariance, W and depend only on coordinate differences. We shall 
take them to be x,y,X — Y, where we used the definitions of Eq. (2.9) and assumed analogous 
definitions for the yi,y 2 variables. Thus 

G^^KxuX 2 ,yi,y 2 ) = G^^\x,y,X- Y ), (2.20) 

W{xi,X 2 ,yi,y 2 ) = W{x, y,X-Y), (2.21) 

where we retained the same symbols to keep the notation simple. Equation (2.18) becomes 

G(4)(x, yX-Y) = 8 i[X-YYy 2 {x- y)] 82 [X-Y-y^{x- y)] 

+ J d^zd'^Zd^z' d^Z' 81 [X - Z + r] 2 {x - z)] 82 [X - Z - yi{x - z)] 

X Wiz, z'; Z - Z')G^'^^ {z\ y; Z' - Y) (2.22) 

We define the Fourier transforms of G, W, 81 and 82 , 

= J d'^x d'^y d‘^{X — Y) G^'^\x,y, X — Y) exp{ipx — ip'y) exp [iQ(X — F)] , (2.23) 
W{p,p';Q) = J d'^xd'^yd'^{X — Y)W{x,y;X — Y) exp{ipx — ip'y) exp [iQ(X — F)] , (2.24) 

and 

S^{p) = J d^z 8 j{z) , (2.25) 

with 8 j (p) being the momentum-space propagator for Xj ■ From the above, we deduce the relation 
between the conjugate momenta of xi,X 2 , which we shall call herepi,p 2 ) and the conjugate momenta 
of X, X, denoted above as p, Q: 

Q=Pi+P 2 , P = V 2 Pi-ViP 2 , (2.26) 

Pi=iliQ+P, P 2 = r] 2 Q-P- (2.27) 
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Analogous relations hold between the conjugate momenta of j/i, 2/2 and those of y, Y. 

For convenience, we also define 

S{p\Q) = Si{'qiQ+p) S2{p2Q-p) ■ (2.28) 

We may now rewrite the Dyson-Schwinger Eq. (2.22) for the 4-point function, in momentum space 

&^\p.p'-.Q) = {2n)^5\p-p')S{p-Q) + S{p-Q) J ^ W(p,fc;Q)GW(fc,p';Q). (2.29) 

We shall use Eq. (2.29) to derive the Bethe-Salpeter equation for the wavefunctions of Sec. 2.1. 


2.3 Completeness relation and decomposition of the 4-point function 

To compute the Bethe-Salpeter wavefunctions of Sec. 2.1, we have to decompose the 4-point Green’s 
function of Sec. 2.2 using the one- and two-particle completeness relation. Then, Eq. (2.29) will 
yield the equations which the wavefunctions and 4>Q^q satisfy. 

Including the one- and two-particle states with the same quantum numbers as xi and X 2 , the 
completeness relation is 



<PQ 

(27r)3 2a;Q,„ 


\Bci,u){Bci.n\ + 


(fq (PQ 1 

(27r)3 (27r)3 2a;Q,q 2eQ,q 


|WQ,q)(WQ,ql, 


(2.30) 


where we have assumed the standard relativistic normalization of one-particle momentum eigen¬ 
states (p|k) = 2£'p (27r)^(5^(p — k), with Ep being the energy of the state |p). To lowest (zeroth) 
order in the interaction strength. 


2wQ,q2eQ,q ~ 2Ei(q; Q) 2i?2(q; Q) , 


(2.31) 


where 

^^i(q;Q) = + + ^1 , A 2 (q; Q) = -y/ {mQ, - q)^ + ml . (2.32) 

Next, we insert the unity operator of Eq. (2.30) in to obtain the decomposition 

G(4) (x, y; a - F) = ^ GW (cr, y; A - F) + G^/^ (x, y,X-Y), (2.33) 

n 


where Gn'^ (x, y; A — F) and G^^ (x, y; A — F) are the contributions of the bound and the scattering 
states, respectively. We compute them below. We shall make use of the fact that a non-zero 
contribution to G^"^^ from a one- or two-particle state arises only when two annihilation operators 
act on that state to obtain the quantum numbers of the vacuum. Moreover, in order to extract the 
poles and the branch-cuts of G^"^^ and G ^^, we will use the integral representation of the 0-function. 


0 {z) 


i 

27r 



^—ikz 

k + ie ^ 


(2.34) 


and 

e [min(x?,X 2 ) - max(y5’,y^)] = 0 [X° - Y° + /i_(x°) - h+{y°)] , (2.35) 

where 

h±{x°) = ^{r]2-r]i)x° ±^\x°\ . (2.36) 

(For Eqs. (2.35) and (2.36), see Ref. [68] and appendix B.) 
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Contribution of the bound states to the 4-point function 

The contribution of the nth bound state to is 


G^:t\x,v,X-Y) = 
f (fK 1 


(27r)3 2a;K,r 

(PK 1 

(27r)3 2a;K,r 

(Pk 1 


{n\Txi{xi)x2{x2)\BK,n){BK^n\Txl{yi)xl{y2)\^) d [min(a;?,a:2) - max{y°,y2)] 

e [XO - yo + h_{x°) - h+{y°)] 


(27r)3 2 wk,,i 

.0 exp {-i [K° - ccK,n] [X° - y° + - h+{y°)] } 


27r 


dK^ 


— oo 
74 




K° - iOK^n + ii 

exp {-i [K° - WK.n] [h_{x°) - ft.+ (2/°)]} 


(27r)4 “ 2a;K,n [JfO - WK,™ + ie] 

where in the third step we made use of the integral representation of the 0 function, given in 
Eq. (2.34), which introduces the integration over K^. The Fourier transform of the above is 


&nHp,p;Q) = J d*xd%d\X - Y) Gl^\x,y,X- Y) 


= ^^Ax)^hAy) 


exp {-i [Q° - u;Q,n] [h-{x°) - h+{y°)]} 
^^Q,n ^Q,n 4” 


At <5° —>■ (jJQ^n, this becomes 


2wq_„ — wq_„ + ie\ 


. (2.37) 


(2.38) 


Equation (2.37) is the contribution of the nth bound state to &P{p,p'; Q). Evidently, the scattering 
amplitude has a pole at energy equal to the bound-state energy. 


Contribution of two-particle scattering states to the 4-point function 

Following similar steps, we find the contribution of the two-particle states to G*-^\ 


G\^\x,y-X-Y) = 
f d^k d^K 


(27r)3 (27r)"^ 


e-^Kix-Y) 


exp {-i [K° - cjK.k] [fe-(a;°) - h+{y°)]} 
2a;K,k 2eK,k [K^ — WK.k + *e] 


The Fourier transform of G^\x, y\X — Y) is 


g\j\p.p'\Q) = 

■ f d?q 


(27r)3 


d^xd% $Q,q(x) $^,q(y) 


exp {-i [Q° - WQ^q] [h_{x°) - h+{y°)]} 


2 WQ,q 2 £q q [Q^ ^Q.q 4~ 


(2.39) 


Clearly, the contribution of the two-particle states to G*^^^(p,p'; Q) gives rise to a branch-cut in the 
scattering amplitude. 

Summing the contributions from the bound and the scattering states, we obtain the decompo¬ 
sition of the momentum-space 4-point function 

G(4)(p,p'; Q)=Y^ GG){p,p'- Q) + &^\p,p'-Q) . (2.40) 
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We shall now combine the Dyson-Schwinger Eq. (2.29) and Eq. (2.40), to obtain the Bethe-Salpeter 
equation for the wavefunctions „ and $Q.q. 


2.4 The Bethe-Salpeter equation for bound and scattering states 

We introduce the operator 


A{p,p'-,Q) 


{2tt)^S‘^{p - p') 

S{P-.Q) 


W(p,p'-Q). 


Then, the Dyson-Schwinger Eq. (2.29) can be cast in the form 

/ (^7T^ MP,k;Q) &^'>ik,p;Q) = (2TTfS'^{p - p) . 


This is formally solved by 

&^\p,p-, Q) = E ^ Q)Clip'-, Q) 

Cn(W) 


Fa{p\ Q)Fl{p-, Q) , 


(2.41) 


(2.42) 


(2.43) 


where Cn{p\ Q) and Fa(p\ Q) are the eigenfunctions of the discrete and the continuous spectrum of 
the operator A{p,q;Q), with eigenvalues c„(Q) and faiQ) respectively. 


f d^k 

J ( 27 r )4 

f d^k 

J ( 2^)4 


A{p, k] Q)Cn{k; Q) = Cn{Q) Cn{p] Q) , 
A{p, k; Q)Fa{k; Q) = fa{Q) Fa{p; Q) , 


(2.44) 

(2.45) 


normalised according to 


E Cn{p\ Q)Ci{p'- Q)+J da Faip; Q)Fl{p'-Q) = (2^)" 5\p - p') . (2.46) 


We may now collect Eqs. (2.37), (2.39), (2.40), and (2.43). Matching the various 
between (2.40) and (2.43), we deduce the following. For the discrete spectrum: 

contributions 

Cn{p; Q) (X J 

d^x 4'q.„(x) 

] h-{x°) 

5 

(2.47) 

cl(p'',Q) cx J 



(2.48) 

Cn{Q) oc 1 

— ^Q,n/Q^ ■ 


(2.49) 

For the continuous spectrum, we identify a — >■ q, and deduce 



Fa{p;Q) oc J 

d*x 4>Q,q(a:) 


(2.50) 

f1{p';Q) oc j 

A ^Q,q(a:) 

.]h+iv°) 

•> 

(2.51) 

faiQ) oc 1 

— a;Q^q/Q° • 


(2.52) 


The relations (2.47), (2.48), (2.50) and (2.51) are stipulated because c„,/o are independent of the 
momenta all the p,p'-dependent factors must arise from the eigenfunctions, C„ and Fa- The 
relations (2.49), (2.52) are warranted so that Cn and Fa are not singular in the limit —)■ 

and —>■ wq.q respectively; the factors [1 — [1 — WQ^q/Q°]“^ cannot be part of the 

eigenfunctions, and thus belong to the eigenvalues. 
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Inserting the above into the eigenvalue equations (2.44) and (2.45), and taking the limits —)■ 

and —>■ wq.q respectively, we obtain the Bethe-Salpeter equations for the bound and the 

scattering states^ 

^Q.n(p) = S'(p;Q) j W{p,k;Q) ^Q^nik) , (2.53) 

^Q,q(p) = Sip; Q) / Wip, k; Q) $Q,q(fc) . (2.54) 

These are homogeneous equations and do not determine the normalisation of „ and 4>q.q. 
Moreover, because we do not know the exact eigenvalues c„ and fa, we cannot use Eq. (2.46) to 
obtain the normalisation of 'I'q.n and 4>q.q. We derive their normalisation in the next section. 


2.5 Normalization of the Bethe-Salpeter wavefunctions 

We derive the normalisation of the wavefunctions Tq „ and 4>q_q from the inhomogeneous Dyson- 
Schwinger Eq. (2.29), or equivalently from Eq. (2.42), using the method described in Ref. [68]. 

We define the symbolic multiplication 

d^k 


[Or 02]ip,p';Q) = 


(27r) 


Oiip,k;Q) 02ik,p';Q ), 


(2.55) 


and the unity operator I{p,p') = (271)“^ (5'^(p — p'). Then, Eq. (2.42) can be expressed in symbolic 
form 

= &^'^ A = I. (2.56) 

We differentiate Eq. (2.56) over and re-use it, to obtain 


~(4) dA -(4) ^ d&‘^'^ 

dQO dQO 


(2.57) 


We shall use Eq. (2.57) to obtain the normalisation of the Bethe-Salpeter wavefunctions. For later 
convenience, we define 

'dAip,p';Qy 


Nnip,p';Q) = i 

Nq{p,p;Q) = i 


dQ° 

dAip,p';Q) 
dQ° 


- Q^—LOQ,rL 

- QO— LjQ^q 


and their Fourier transforms, 
Nn{x,x';Cl) = i 


Nq(x,a:';Q) = i 


d 

r d'^p d^p' 

Aip,p';Q) e^P'^' 

[dQ°j 

(271)"* (27r)4 

d 

f d'^p d'^p' 

e-^P^ Aip,p';Q) 

[dQ° J 

{2TTy (27r)^ 






(2.58) 

(2.59) 

(2.60) 
(2.61) 


Bound states 


Substituting the contribution to the 4-point function from the nth bound state, given in Eq. (2.38), 
into Eq. (2.57), and taking the limit —>■ ujq^u, we obtain the normalisation condition 

f d^T) d^v^ ^ ^ ~ 

J (27r)'^ ^27 i-)4 ^nip,P ;Q) 'I'q_„/(p ) = 2a;q_„ . (2.62) 

In coordinate space, this becomes 

y d^a;d‘‘a;'4'Q_„(a;) N„(a;,a;';Q)4'q^„/(a;') = 2wq_„ (5„„/ . (2.63) 

^It is possible to obtain Eq. (2.53) more easily, by taking the residue of both sides of Eq. (2.40) at ^Q,n- 

However, this is not possible for the two-particle states. 
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Two-particle states 

Substituting the contribution to the 4-point function from the two-particle states, Eq. (2.39), into 
Eq. (2.57), we deduce the normalisation condition^ 

In coordinate space, this becomes 


J d^x d*x' $Q,q(a;) Nq{x, x'] Q) $Q^q/ (x') = 2wQ,q 2eQ^q (27r)3(53(q - q') . 


(2.65) 


Note that in the fully relativistic case, the normalisation of the wavefunctions depends in general 
on the potential. 


2.6 The instantaneous approximation and the Schrddinger equation 

In the non-relativistic regime, it is possible to simplify the Bethe-Salpeter equations. The momen¬ 
tum exchange between the two unbound particles is |q| ~ /rurei, while between two bound particles 
|q| ~ /io, with a characterising the interaction strength; in either case, for a,Urei «C 1, the energy 
exchange is ~ q^/(2p,) ^ |q|. It is then reasonable to ignore the dependence of the kernel 
W(p,p'-,Q) on This is the instantaneous approximation.® In fact, in the cases of interest, 

W(p,p'-,Q) depends only on |p — p'|, rather than on p and p' separately, and it does not depend 
on Q (except perhaps for Q^, which, in the non-relativistic regime, we shall approximate with 
~ TO^). We shall thus assume that 

W(p,p';Q)~ W(Ip-p'I). (2.66) 


In this approximation, we deduce from the Bethe-Salpeter Eqs. (2.53), (2.54), that 'I'Q^„(p)/S'(p; Q) 
and $Q^q(p)/S'(p; Q) are independent ofp°. We define 


5o(p;Q) 

^Q.n(p) 

0Q,q(P) 




So{v\Q) ' 
S{p-,Q) . 

' 5o(p;Q) 

2^Q,q ['5 '(p;(3) 


y^2A/Q(p) 

/2W'q(p) 


^Q.q(p) > 


where we choose the normalization factor 


A^q(p) = 


^i(p;Q)^2(p;Q) 
■Ei(p;Q) + E2 (p;Q) ’ 


(2.67) 

( 2 . 68 ) 

(2.69) 


(2.70) 


such that we recover the conventional normalisation for V'Q.n and ((>Q.q, as we shall see in Sec. 2.7. 
We calculate 5o(p;(5) in appendix C. Multiplying both sides of Eqs. (2.68), (2.69) with S{p-,Q), 
integrating over and using Eq. (2.67), it follows that 

^Q,n(p) = ^j2AfQ{p) J ^ ^Q.n(p) = ^J2^fQ{p) J Sx 4'q,„({x° = 0, x}) , (2.71) 


<^Q,q(p) 


^A/q (p) 

2eQ,q 


dp^ 

27r 


®Q,q(p) 




X $Q,q({a; 


0 


0,x}) 


(2.72) 


“^In fact, from Eqs. (2.39) and (2.57) we obtain the normalisation condition described in Eq. (2.64) with the functions 
^Q.q(p) replaced by Eq(p;Q) = f ^Q,q(3:) (c.f. Eq. (2.50)). Then, taking ^Q.qi 

we obtain the exact form of Eq. (2.64). 

^For a discussion on relativistic corrections, see e.g. Ref. [69] and references within. 
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V'Q,Tt(p) and </'Q,q(p) are sometimes called the “equal-time wavefunctions”. From the above, and 
recalling Eqs. (2.7), (2.8) we see that 

^Q,n(p) = '0Q,n(P) . ^Q,q(P) = ^Q,q(P) • (2-73) 


Given the definitions (2.68), (2.69) and Eqs. (2.71), (2.72), the Bethe-Salpeter Eqs. (2.53) and 
(2.54) become 


^Q,n(p) ^ 

a/2A/'q(p) 5o(p;Q) 

^Q,q(p) ^ 

a/2A/'q(p) So{p;Q) 

Non-relativistic approximation 


d'^k yV(|p —k|) 
(27r)3 ^2A/'Q(k) 

d^k >V(|p-k|) 
(27r)3 y'ZSfdik) 


^Q.n(k) , 

with Q° = , 

(2.74) 

^Q,q(k) , 

with Q° = WQ,q . 

(2.75) 


Using the non-relativistic approximation described in appendix C, Eqs. (C.18) - (C.22), and setting, 
in accordance to Eq. (C.21), 


ujQ^n = m + Q'^/2m + £n , (2.76) 

WQ,q = TO-I-Q^/2m-I-fq , (2.77) 

equations (2.74) and (2.75) become 

(-|; + Mp) = WdP - kl) <»4k) , (2.78) 

(-1; + £,) 0,(P) = - ^ / W(|P - k|) 0,(k). (2.79) 


These are the Schrddinger equations for the bound and the scattering states in momentum space. 
They are eigenvalue equations, and as such, their solutions determine £n and £q- Because in 
Eqs. (2.78), (2.79), all dependence on the CM momentum Q has been eliminated, we have dropped 
this subscript from the ijj, (f> wavefunctions, but kept the same symbols in order to avoid cluttering 
the notation. Note that from Eq. (2.76), it follows that the mass of the bound state is 


Mji — TO -f £n • 


(2.80) 


It is convenient to Fourier-transform Eqs. (2.78) and (2.79) to coordinate space. We set 

d^p 


^n(l’) = 

= 


(27r)3 
• d^p 
{2 ttY 


Ynip) , Ynip) = J d^r V’n(i’) e 

<^q(p) e'P '’, (^q(p) = [ dYr Yqir) 


(2.81) 

(2.82) 


Acting on both sides of Eqs. (2.78), (2.79) with f we obtain the Schrddinger equations 

in coordinate space 


- 2 ;;+ 

■ 


1pn(r) = £nYn(r) , 
'/^q(p) “ ^q^q(p) j 


where V (r) is the non-relativistic potential, 


U(r) = - 


1 r d^k 
iAmfi J (27r)3 


W(k) 


^ik-] 


(2.83) 

(2.84) 

(2.85) 


We quote the bound-state and scattering-state solutions of the Schrddinger equation for a Coulomb 
potential, in appendix F, and use them in our computations in Sec. 5. 
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2.7 Normalization of the Schrodinger wavefunctions 

To find the normalization of tpn and d)p, it is easiest to follow a similar procedure to that of Sec. 2.5. 
We first define 


GuHp,p';Q) 

GuHp,p';Q) 


dp° dp'° , 


'dp" dp'"^(4) 


27r 27r 




( 2 . 86 ) 

1 »^n(p)C(p') ^2 87) 

y/2A/q(p)2A/q(p') ^ujQ^n [Q° — WQ,n + *e] 

_1_ f d^q »<^q(p)^q(p') 

\/2A/q(p)2A/q(p') J (27r)3 2a;Q^q + ie] 

( 2 . 88 ) 


Integrating Eq. (2.29) with respect top°,fc° yields 

^^'‘^(P,P';Q) = (27r)^d^(p- p')5o(p;Q) +5o(p;Q) J >V(|p - k|) 0(^)(k,p';Q). (2.89) 


Following the steps of Sec. 2.5, we find the equivalent of Eq. (2.57), 


r d^k 

J (^ 


0^"^(p,k;Q) 


d 


5o(k; Q) 


g^^HKp'-,Q) 


-^^(^Hp,p';Q). 


From Eq. (C.14), 


d 


.5o(k;Q)_ 


-i2AfQ(k) 2Q". 


We may now obtain the normalisation conditions for the wavefunctions and ^q. 


(2.90) 


(2.91) 


Bound states 

Close to the pole, at —?> wq^„, we may substitute the contribution from the nth bound state, 

Eq. (2.87), into (2.90). We obtain 

/ ^ ^ ^ 

Two-particle states 

Substituting the contribution from the two-particle states, Eq. (2.88) into (2.90), and for —>• wg^q, 

we deduce the normalisation condition 

/ ^ ^ (27r)^d3(q - q') ■ (2-93) 


Note that in obtaining the normalisation conditions (2.92) and (2.93), we did not make use of 
the non-relativistic expansions of the factors A/g(p) and £g,q, given in Eqs. (3.32) and (3.34). 

3 Radiative level transitions 

In this section, we determine the amplitudes for the radiative BSE and de-excitation processes 

Xl + X2 —>■ (XlX2)bound + <P ) (3.1) 

(XlX2)bound,n' ^ (XlX2)bound,n + ‘p t (3-2) 
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rjiK + k 


rj2K - k 




riiP + p 


G(5) / 


mP -p 


rjiK + k 


t] 2K — k 


(5(4) 




ViP + P 




\ ■^2^’ - P 

^(5) 


I 


Figure 3. Diagrammatic representation of the equation (3.12) for the 5-point function (5*-®). stands 

for the full propagator of the force mediator ip, which may be either a scalar or a vector boson. 


in terms of the bound-state and scattering-state wavefunctions computed in Sec. 2 and a pertur¬ 
bative interaction which describes the emission of the force mediator. The S-matrix elements of 
interest are 


out(^P,n i |^K.k)in — (^P,n : I ^ I ^K,k) ; (^■^) 

out(Sp.n ; Pp^ I B'K,n')in = {Bp ,n Pp^ \ S | , (3.4) 

where the indices stand for the momenta and the quantum numbers of the corresponding states, as 
stated in the beginning of Sec. 2. 

3.1 The 5-point Green’s function 

Since Bq^u and f^q.q are generated by the action of Xi and xl vacuum (c.f. Eqs. (2.4), 

(2.6)), in order to compute the S-matrix elements of Eqs. (3.3) and (3.4) we need to consider the 
5-point function 

G^^\x^,xi,X2 ; yi,y2) = {fl\T(p{X^)xiixi)x2ix2)xi{yi)xUy2)\^) ■ (3.5) 

We define the Fourier transform 

G('’^(P^,Pi,P2; ki,k 2 ) = 

= j d^X^d‘^xid^X2d^yid% G^^\x^,xi,X2; yi,y2) ■ (3.6) 

As in Eq. (2.9), we set 


X = xi — X2 , X = riixi + r]2X2 (3.7) 

y = yi-y 2 , Y = yiyi + 772^2 , (3.8) 


and rewrite the above as 

{P^, r]iP + p, r]2P - P; ipK + k, 772 AT - k) = 

[ d'^X^d'^Xd'^xd'^Yd'^y e^iP^^^.+PX+px-Ky-ky) + y 2 X,X-yix; Y+ r] 2 y,Y-r]iy), 

(3.9) 
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i.e. the conjugate momenta of X, x are P^p, and the conjugate momenta of Y, y are K, k defined as 


P = pi+P2 , p = r]2Pi - riiP2 , (3.10) 

K = ki + k2 , k = r]2ki - r]ik2 ■ (3-11) 

The 5-point Green’s function ,xi^X 2 ]yi,y 2 ) is equal to the sum of all connected dia¬ 

grams with five external points. The momentum-space is sketched in Fig. 3, and can be written 
as 


/ ^4 / 

(27r)4(54(x -P-P^) mP + p', mP - p'; mK + k', y2K - k') &^\k\ k; K) 

where 


is the (fi propagator, with 


q (p \ _ iZ^(Pcp) 

P^-ml+ze 


Z^(q) = |(n|(^(0)|v?q)p 


being the field-strength renormalisation parameter for p. is defined via the relation 
iC^^\P^,Pi,P2; ki,k2) = S^{P^)Si{pi)S2{p2)Siiki)S2{k2) iA^^\P^,pi,p2] ki,k2), 


(3.12) 

(3.13) 

(3.14) 


(3.15) 


where 


iPijP 2 ] ki,k 2 ) = sum of all connected diagrams. 


(3.16) 


Note that may include diagrams that are not fully connected, i.e. diagrams in which external 
legs are disconnected from each other,® but it does not, of course, include vacuum bubble diagrams. 
(If only fully connected diagrams contributed to then A^^'^ would simply be the sum of all 
connected and amputated diagrams, as conventionally defined.) For later convenience, we also 
define as the sum of all connected diagrams with only the (/3-leg amputated. 


iC^^\P^,Pi,P 2 ; ki,k 2 ) = S^{P^) 


,{P<p,Pi,P 2 ; ki,k 2 ) ■ 


Then, A^^'^ appearing in Eq. (3.12), becomes 


(3.17) 


iA^^\Py,, riiP+p,r]2P-p] r]iK + k,r]2K-k) 


iC!p-s,mp(Pv’'niP + P^mP - p; ViK + k,y2K -k) 
S{p-,P)S{k-,K) 


where we remind that S{p;P) = Si{r]iP + p) S 2 iri 2 P — p) (c.f. Eq. (2.28)). We sketch Eqs. (3.15) 
and (3.17) in Fig. 4. 


3.2 Transition amplitudes 

We now extract the S-matrix elements of Eqs. (3.3) and (3.4) from the 5-point Green’s function of 
Eq. (3.5). Our analysis follows closely Sec. 7.2 of Ref. [70]. 

Let us first focus on the BSF amplitude of Eq. (3.3), for which 

P° ^ W^(P^), P° ^ C^P.n, K° ^ WK.k . (3.19) 

®In fact, the lowest-order contribution to for the transition processes considered in Sec. 5, arises from diagrams 
that are not fully connected, as shown in Fig. 7. However, the entire transition processes are described by fully 
connected diagrams, shown in Fig. 8. 
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Figure 4. Diagrammatic representation of the equations (3.15) and (3.17). stands for the sum of all 
connected diagrams with no legs amputated; this includes not fully connected diagrams. is equal 

to with only the ip-leg amputated. —#— and 0 stand for the xi s-nd X 2 full propagators, 
respectively. stands for the full propagator of the force mediator ip, which may be either a scalar 

or a vector boson. 


In this limit, the Lehmann-Symanzik-Zimmermann reduction formula yields 


J J d^X J d^Y G^^\X^,X + r,2 X,X- rjix; Y + my^Y - r,iy) ~ 


f (ll|V3(0)|v?p^) 

i {^\Txi{'n2x)X2{-pix)\Bp^ny 

2 uj^{V^){P^ - w,^(P^) -b ie) 

2 cup_„(po — cup^ra-b ie) 


d^k' i{ltK,k>\Tx\{v2y)xl{-Viy)\^) 


(27r)32e: 


K,k' 


2 a;K,k' — WK.k' + te) 




K,k' 


(3.20) 


Here, the ^ sign means that the two sides have the same singularities in the limit (3.19); to 
compute the S-matrix element, we need to extract the residues of these singularities from both 
sides of Eq. (3.20). 

In the above expression, the correlation functions involving the xi X 2 fields correspond to 
the bound and scattering state wavefunctions (c.f. Eqs. (A. 5), (A. 8 )). The correlation function 
involving the ip field is the p field-strength renormalisation parameter (c.f. Eq. (3.14)). We Fourier- 
transform Eq. (3.20) with respect to x,y, to obtain 




i ./zyjpy) 

i^-p,n{p) 

2 uj,p{P,p){P° - -b ie) 

2 u)p^n{P^ — Wp,n + ie) 


d^k' 


*^K,k'(9) 


(27r)32eK,k' 2ajK,k'(Ar° — cuK.k' + *e) 


— (Sp,„,7’pJS|Z^K.k'). (3.21) 


The left side of the above equation is [P^, yiP +p, r] 2 P — p; piK -b q, 772 AT — q), which may be 
decomposed according to Eq. (3.12). Recalling Eqs. (2.37), (2.39) and (2.40) for the 4-point function 
and keeping only the leading singularities in the limit (3.19), the left side of Eq. (3.21) becomes 


{P,f, piP + p, r] 2 P -P', t]iK + q, q 2 K - 9 ) ~ 

r d^p' d\' i^P,n{p)^*p^n(p') f d^k' ^^K,k'(g')^K.k'(g) 

P^-rrP^ + ie J (27r)4 (27r)4 2 wp,„ [P° - wp,„ -b ie] J (27r)3 2 wK,k' 2 £K,k'[AT^ - WK.k' + ie] 

X { 27 t)^S\K -P-P^) iA^^\P^, piP + p', P2P - P\ PiK + q', P2K - q') . 
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At —)■ oJipiV^) and > a;p,„, this expression has the same poles as the right side of Eq. (3.21). 

Identifying their residues, we obtain 



( , f d^k' »^K.k'(g')^^c.k'(g) 

J (271)"* (27r)4 P,n P J (27r)3 2a;K,k'2eK,k'[Ar° — cjK.k'+ *e] 

{2tt)‘^S'^{K - P - P^) iA^^\P^,riiP + p ,r] 2 P - p-, piK + q ,r] 2 K - q) 


(Pk' 


*®K,k'(9) 


(27r)32eK.k' 2wK,k'[Ar° — WK.k' + ie] 




K.k' 


We still have to extract the leading singularity at t^^K.k- We multiply both sides of the 

above expression with A^k(9, <7^^ K)$K,k(9"), integrate over q and g", and use the orthonormality 
condition (2.64), to obtain the S-matrix element for BSF 


(Sp,„;^,JS|7^K,k) = y^^(P^) I ^K.k(g) 

X {2Tr)^S\K -P-P^) {P^, ViP + P, mP - p; mK + q, q 2 K - q) . (3.22) 

Following similar steps, we obtain the S-matrix element for transition between discrete energy levels, 

{Bp,n ; I S I ^K.n')- 

In standard notation, we write the S-matrix elements as 

(ep,„ ; ¥>p J S I Z^K.k) = iAirf 5\K-P- P^) iMk^n , (3.23) 

(Sp,„ ; J S I ^K.n') = (2^)" S\K -P-P^) , (3.24) 

with 

= ^Z^(P^) / ^P.n(p) ®K,k(<z) 

X A^^^ {P^, r]iP + p, r] 2 P - p; pi AT -|- q, 772 AT - q) , (3.25) 


= ^Z^(P^) / ^P.n(A) ^'Kx(<z) 

X {P^, r]iP + p, 772 P - p; 771 AT -I- q, 772 AT - q). (3.26) 

If non-fully connected diagrams contribute to the perturbative part of the transition amplitudes, 
then should be replaced by using Eq. (3.18). In this case, we obtain 


■Mlc^n = 



d^p d'^q ^p,„(p) $K.k(9) 

(27r)4(2^)4 S{p-,P) S{q;K) 

c!p-a.mpiPvAiP + p,mP - p; mK + q,mK -q), (3.27) 


M 



d^p d'^q ^'p,„(p) ^K.n'(9) 

(2^)4 (2^)4 S{p-P) S{q-K) 

^C^^-arnp{PvAiP+ PA 2 P-P', mK + q,q 2 K - q) . (3.28) 


In the case of a vector mediator Zy, becomes the charge-renormalisation parameter and the 
amplitudes contain the polarisation vector e^, i.e. A^^'^ = and C^^2amp = 
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3.3 Instantaneous approximation 


In the instantaneous and non-relativistic approximations, we may express the transition amplitudes 
in terms of the Schrodinger wavefunctions defined in Eqs. (2.68), (2.69), as follows 




^ - f (Pp (Pq <^k(q) 

^ ^ (27r)3 (27r)3 ^27Vp(p) 2AfK(q) 


-^trans(qjP) ; 


r (Pp (pq ^n(p) ^n'(q) 

J (27r)3 (27r)3 ^2Afp(p) 2AfK(q) 


.^^trans (q? P) ^ 


(3.29) 

(3.30) 


where we took Z^(P^) ~ 1 to lowest order, and set 

= go(q;g)5„(p;P) /-ftP-p; »./?+<;■ 


(3.31) 


It is sufficient for our purposes, and consistent with our approximation (see footnote 14), to 
expand the normalisation factors up to first order in p^, q^, as follows 


A^q(p) - M 


\ + A 


2 p? 

V "lyj 


a/ 2A/p (p) 2 A/k (q) 2/r 


r p' + q" 

1 

CO 

re 


V ^ y 


(3.32) 

(3.33) 


The p^, terms in Eq. (3.33) introduce corrections of order P and (see appendix E), where 
a parametrises the strength of the interaction (for a Coulomb potential, it is the fine-structure 
constant) and gives the expectation value of the relative velocity inside the bound state. Similar 
corrections arise also in Altrans (see appendix E). We shall retain such corrections only where the 
dominant term in the respective expansion vanishes, as in the case of degenerate particles interacting 
via scalar boson exchange (see Sec. 5.1). Moreover, from Eqs. (2.31) and (2.77), we find that to 
zeroth order in the relative velocity, 

EK.k —M- (3.34) 


Because £K.k factors out of the integrals, as seen in Eq. (3.29), we neglect corrections, which 
always produce subdominant terms in In addition, in our computations, we consistently ignore 
corrections that involve at least one power of the total momentum of any of the initial or final states 
(denoted typically with capital letters). In the CM frame, these momenta are of order ~ 
which renders their scalar products with any other momenta, of higher order in a and t;rei than the 
p^,q^ corrections. 

We shall employ Eqs. (3.29) - (3.33) to evaluate the transition amplitudes of Sec. 5. 


3.4 On-shell approximation 

Let us now consider the case when the perturbative part of the transition amplitude consists 
only of fully connected diagrams.^ Then, A^^\Pip,r]iP+p, rj 2 P — p', ViK + q, rj 2 K — q) is the pertur¬ 
bative amplitude for the 2 —>■ 3 transition (with no on-shell conditions imposed). Equation (3.31) 
becomes 


.^trans(q5 P) 


fdp°dq° Sip-P)S{q-K) 

J 27r 2 tt 5o(p;P)5o(q;if) 


A^^\P^,r]iP +p,ri2P 


p\ r]iK + q,p 2 K - q)- (3.35) 


Provided that A^^'^ has no singularities in and g®,® the integrations over p^,q^ force the 
evaluation of on the poles of S{p-P) and S{q;K) that are located in either the lower or 

^This is the case if the initial-state particles and the particles participating in the bound state are different. 

p,7]2P — P'i — O') may have singularities in p^ and (for the energy of interest, 

= cjK.k)? if fh® initial-state particles and the particles participating in the bound state are different. 
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upper and complex planes (depending on the choice of integration contours). As described 
in appendix C, each integration picks out two poles: One physical pole, which corresponds to 
setting one of the particles on-shell, and one unphysical pole, where the energy of the other particle 
is negative (c.f. Eqs. (C.ll), (C.12)). In the non-relativistic regime, the contribution from the 
physical pole dominates. For concreteness, let us take these poles to be in the lower 9 ° complex 
planes (as in Eq. (C.12)), 

pO =+ Ei{p;P) - le, 

= -r^,K^ + E,{q;K)-ie. 

Fixing and to the pole values means that the energies of the Xi,X 2 particles in the bound 
state and in the two-particle states are specified as functions of the 3-momenta p,q, P,K and the 
quantum numbers n and k (note that = uip^n and = wk.Ic), as follows 

p? = -f = E;i(p;P) , = a;p,„ - Ei(p;P) ~ £; 2 (p;P)- f - p^/(2/r), (3.36) 

= + = Ei(q;K), = cuK.k - Ei(q; K) ~ E 2 (q; K) + fk - qV(2M), (3-37) 


where we used Eqs. (2.76), (2.77) and (C.18), (C.19). Evidently, in both the bound state and the 
two-particle state, the xi degree of freedom is off-shell, by £„ — p^/ 2 p and E\^ — q^/ 2 p respectively. 
However, (p^)/(2p) ~ and (q^)/(2p) ^ provided that En,E^ ^ p,® we may ignore this 
small deviation from the on-shell condition and evaluate on-shell. Then, from Eq. (3.31), we 
obtain 


Altrans(Qi P) 


+ P,V 2 P 


p; riiK + q,r]2K - 


q) 


on—shell 


(3.38) 


This is the approximation presented in Sec. 5.3 of Ref. [70]. Note that, for consistency, when 
using Eq. (3.38) inside Eqs. (3.29) and (3.30), the normalisation factor of Eq. (3.33) should be 
approximated to zeroth order in p®, q^. Indeed, corrections of the order p®, q^ arise not only due to 
the normalisation factor, but also due to the off-shellness of the amplitude When important, 

such corrections should be included self-consistently, by making use both of the full expansion of 
Eq. (3.33), and of the off-shell momenta of Eqs. (3.36), (3.37) instead of the on-shell conditions. 

We will not make use of Eq. (3.38) in our computations in Sec. 5. 


3.5 Bound-state formation cross-sections 

In the CM frame (K = 0), the differential cross-section for radiative BSE, UK=o,k —t Bp,n + P_p, is 

1 |P| 


dEl 


2 ^ — m\ — ~ 4 m®m| 167r® y/s 


\M 


k—>-n| 5 


(3.39) 


where s = k — + fk)^ (c.f. Eq. (2.77)), and Alk-j-n is found from Eq. (3.29). The 

bound-state and mediator momenta are 


IPI = 


{s-Ml- mlf - 


As 


1/2 


~ {Ek - Er, 


1 - 


mt 


(£k - En 


1/2 


where we used M„ = m + En (c.f. Eq. (2.80)). In addition. 


(s — m\ — m^)'^ — Arnfm^ — Amk = AmfiVrei ■ 


Then 


(£k - fn) 


dEl GAn'^rrPpVrei 


1 - 


m,. 


{Ek - EnY 


1/2 


[A^k^nP ■ 


(3.40) 


(3.41) 


(3.42) 


®For a Coulomb or Yukawa potential with fine structure constant a, these conditions are equivalent to a, <C 1 
(c.f. Eq. (5.4)). 
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Partial-wave decomposition and unitarity 

It will be useful to decompose the amplitude Atk-s-n in partial waves 




2J+1 

47r 


M jPj {cos 0) 


where Pj are the Legendre polynomials, and 

Mj = J dUPj{cos 9) ■ 

Then, Eq. (3.42) gives 


jt") =y^crW 

BSF / y BSF,J 

,7 


with the partial-wave cross-section given by 


J.(n) _ ('^k ^n) 

BSF,J 647r2m^^Urel 


1 - 


m" 


{Sk - Sr^Y 


1/2 


2J+1 

47r 


\M^ 


(3.43) 


(3.44) 

(3.45) 


(3.46) 


Unitarity implies an upper limit on the partial-wave inelastic cross-sections. In the non- 
relativistic regime, for the Jth partial wave, this is [71] 


17inel,J ^ (l^inel.j) u 


(2J-t I)7r 


(3.47) 


For a given inelastic process and associated cross-section, this bound yields an estimate for the 
value of the coupling at which the probability for inelastic scattering saturates. In Sec. 5, we use 
the unitarity bound to deduce the range of validity of our calculation. 


3.6 Bound-state de-excitation rates 

The differential rate for the radiative de-excitation of a bound state, ^p,n + / 2 _p, is 


dr. 


dVL 


327r2M2 


\Mr 


where is found from Eq. (3.30) and jPj is given by Eq. (3.40) with the replacement 

£k —t £n' ■ Setting = m + £n — rn, we obtain 


{£n' - £n) 

i2n'^rrP 


I - 


{SW - Snf 


1/2 


\Mr 


'dn 


(3.48) 


4 Decay of unstable bound states and (co-)annihilation of nnbonnd pairs 

4.1 Non-perturbative amplitude 

If an unbound xi, X 2 pair can (co-)annihilate into a number of light particles fi,---fN^ then 
the xi ~ X 2 bound states are unstable against decay into the same final states (provided that 
this is allowed by angular momentum conservation). For example, unbound and bound particle- 
antiparticle pairs can annihilate and decay, respectively, into force mediators. (Co-)annihilation 
and bound-state decay have the same diagrammatic representation, shown in Fig. 5; the difference 
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Figure 5. Diagrammatic representation of equation (4.5). .. .. represent the full propagators of 

armihilation/decay products fj. 


in evaluating these two processes is which initial state is singled out from the function. In this 
section, we express the (co-)annihilation and the bound-state decay amplitudes, 

out(/l/2 • • •/w I ^K,k)in = (/l/2 ’ ’ ’/at I S | , (4-1) 

out(/l/2 • • •/w I SK,n)in = (/l/2 ’ ’ ’ /at | S | ^K.n) , (4.2) 

in terms of the initial state wavefunction and the perturbative interaction that gives rise to these 
processes. 

To calculate the S-matrix elements (4.1) and (4.2), we need to consider the Green’s function 
Ga^nn{xi,X2 , ...Xn] 2/1,2/2) = {^\T f (xi) f {X 2 ) ■ ■■ f (x n) x{iyi)xUy2)\^) ■ (4.3) 


and its Fourier transform 

W r 

GanniPl, ■■■PN', ^ 1 ,^ 2 ) = II / 

i=l 


d‘^yid'^y 2 e hfciyi-Hfc 2 y 2 ) Gann(a;i,--- ,XN]yi,y 2 ) ■ 

(4.4) 


Let Aannipij''' Pn] ^ 1 ) ^ 2 ) be the sum of all connected and amputated diagrams contributing the 
(co-)annihilation of a yi, X 2 pair with momenta ki,k 2 , into /i,-- - ,/tv particles with momenta 
Pi," ' jPn^^ Then, as sketched in Fig. 5, 


Gann(pi,' ’ ’ ,PN',r]iK + k,p 2 K - k) = 

^ ~ r d'^k' 

= II‘^l(^j) J • • ,Ptv;?7i^ +fc') G('‘^(fc',fc;iF) , (4.5) 

where Sj^ (pj) is the propagator of the fj particle with momentum pj . As always, energy-momentum 
conservation implies that 

Aann(pi, ’ ’' ,Pn', ^ 1 ,^ 2 ) = i{2TT)'^ 5'^ {ki -f “ Pi- Pn) Ml^f{ki,k2-, Pi,-" ,Pn) , (4.6) 

where A4P“/ is the perturbative annihilation amplitude, with no on-shell conditions imposed on 
the incoming and outgoing degrees of freedom. We follow a similar procedure as in Sec. 3.2, and 
determine the S-matrix elements of interest to be 


(/l • • • /tv I S I f^K.k) = {2TTf 5\K-P1---- pn) iMann , (4.7) 

(/i • • • /tv I s I ^K.n) = (2^)^ S^{K - Pi... - P^) iM^l, (4.8) 

already mentioned, for annihilation processes, similar analyses have been carried out in previous works, e.g. [40— 
42, 65, 66, 72-74]. 

^^Note that, in contrast to level transitions processes, all diagrams contributing to the perturbative part of the 
(co-)annihilation processes are fully connected. The amputation of fully connected diagrams is well defined, and 
it is thus sensible to express the full amplitudes for the (co-)annihilation and decay processes of interest, in terms 
of the sum of amputated diagrams. 
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with 


N 


A^ann = n V (P. 

i=i 


d'^q 


— \/2eK.k 


$K,k(g) + q,mK - Ql pi,--- ,pn) 

A^ann (mA + q,r) 2 K -q] Pi, • • • ,pAr) , (4.9) 


I (27r)4 

0k(q) f dq° S{q;K) 


(27r)3 y2A/k(q)y 27r 5o(q;K) 




w 


i=i 


d'^q 

'0n(q) 


^K,n(g) + q,r] 2 K - 9 ; Pi, • • • ,Piv) 

Atann(PlA + q.^K - q] Pi, • • • ,PAr) , (4.10) 


_ ^(g;A) 

(27r)3 i/2A/'K(q) j 271 5o(q;iir) 


where = oiK.k and = ^^K,n, respectively. In the above, Zf-{p) = \{^\fji0)\fj,p)\‘^, with fj^p 
being a fj particle with momentum p. In the second line Eqs. (4.9) and (4.10), we have used the 
instantaneous approximation for the wavefunctions, and set Zf^{p) ~ 1. 


4.2 On-shell approximation 

Following the same arguments as in Sec. 3.3, we may evaluate the perturbative amplitude Alg®))* 
on-shell. This enables us to express the annihilation and decay amplitudes as follows 

Alann 7^ J <^k(q) AIP“‘(77iK -f q, ? 72 K - q; Pi, • • • ,PAr) , (4.11) 

-^™(PiK + q,P2K-q; pi,--- ,pAr) , (4.12) 

where A4 p®(( is the on-shell perturbative (co-)annihilation amplitude. Note that, as discussed below 
Eq. (3.38), the integrands in Eqs. (4.11), (4.12) admit q^ and higher order corrections from the 
normalisation factor of Eq. (3.32) and from the off-shellness of the perturbative amplitude A4P®n . 


4.3 Two-body (co-)annihilation cross-sections and bound-state decay rates 

Let us now focus on the case of decays and (co-)annihilations into two final-state particles. In the 
CM frame (K = 0), the momenta of the final particles are |pi| = IP 2 I = |p|, with |p| ~ WK=o.k/2 = 
(to -|- £k)/2 — to/ 2 in the case of (co-)annihilation, and |p| ~ Al„/2 = (to -I- 5n)/2 — to/ 2 in the 
case of decay. We ignore the masses of the final-state particles for simplicity. The (co-) annihilation 
and decay amplitudes can be expanded in partial waves as follows 

TE / 2 / -I- 1 \ 

A4ann(f2p) = ^ f —- j Pj^{cOS 6»p) Alann/ , (4.13) 

= E Pdcos9p)Mtl,, (4.14) 

e=o ^ 2 

where 

A4ann,^ = ^ dflp (cOS ^p ) At ann (^p) , (4-15) 

J d^2pPr(cos0p)A4':^(llp). (4.16) 
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Here and in the following, the indices in the angle variables specify the vector to which this angle 
refers; a double index denotes the angle between the two vectors. The (co-)annihilation cross-section 
times relative velocity and the decay rate are 


I'^rel — 


^ dec 


fs 


1287r^m/x 

fs f 

647r^TO / 


|A4ann(f^p)Pdf^p = 




1287r2m/i 


E 


21 

47r 


|Xa 






647r^m 


E 


2£- 


47r 


-IMEP: 


(4.17) 

(4.18) 


where fs = lf2 if the final-state particles are identical, or /s = 1 otherwise. We shall now express 
TWann,^ and ^ in terms of the perturbative on-shell annihilation amplitude. 

We expand A4P®a in partial waves, as follows 


‘■ann vH.? 


^=0 


(4.19) 


In general, the expansion coefficients dj may depend on q; in the non-relativistic regime, they can 
be expanded as 

d£(q) ~ Of -f Cyi • q), (4.20) 

where stands for the polarisation vectors of possible final-state vector bosons, and Ti is a 
polynomial function of the scalar products q^ and ea ■ q that vanishes at q = 0. Note that |p| is 
determined by energy conservation, ae and may depend on scalar products such as p^, • p and 

ea' f-B- In the following, we consider only the ai contribution to d^; any corrections arising from 
the q-dependent terms of Eq. (4.20) may be included only in conjunction with similar corrections 
arising from the normalisation factor of Eq. (3.32) and from the off-shellness of the perturbative 
amplitude . 

We may now insert Eq. (4.19) into Eqs. (4.11) and (4.12), and use the formula 


drip Pi! (cos dp 


dfq 

(27r)3 


c(q) |q|^^’£(cosdq,p) = 

, { 2 ^+l)\\ 


iH 2 i+l)e\ 




PlicOsOr) ())k(l 


(4.21) 


and similarly for ■!/)„. We prove Eq. (4.21) in appendix D. Keeping only the q-independent term 
from the expansion of Eq. (4.20), we find 


Ms 


a,\p\^ (2£+l)!! 


M 


(n) 

dec,£ 


{mfiY ^^(2^-I-1)£! dr^ 

g^lpl^ (2^ + 1)!! 

^/ 2 fx{mnY i^{ 2 i-\-\)d. 


r=0 


dilj-Peicos dr)<('k(r) 

' d^ f 

i y 'IIIrEf(cOSdr)'*/'n(r) 


(4.22) 

(4.23) 


Using Eqs. (4.17) and (4.18), we hnd the contribution of the £ partial wave to a’ann'^rei and T^”^ to 
be 


where 


(^ann ^rel)^ — ‘^'^.ann : 

T'(^) C* 

^ £,dec ~ ^i,dec , 

^ [U/(2£)!f fs 

^ 2 £ + 1 327rmfj, ’ 


(4.24) 

(4.25) 


(4.26) 
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and 


Si^ann — 

4^+2 7r2 ^2A 

dP 

q{n) _ 

‘^e,dec — 

[( 2 t + l)\/{£r)^Y 

d^- 

4^+2 7r2 /.t2^ 

dr^ 


dn^Pgicos 0r)^n(r) 


r=0 

2 


r=0 


(4.27) 

(4.28) 


In the limit where the interaction in the two-particle state can be neglected, (/)k(r) = and 
5,.ann = (|k|/M)^^ = r;?il. 

Similar analyses to the above for the non-perturbative annihilation cross-section, have been 
performed in Refs. [65, 66], where also the Sommerfeld enhancement factors of Eq. (4.27) have been 
computed for a Yukawa potential. 


5 Bound-state formation, de-excitation and decay rates for specific in¬ 
teractions 


We now focus on specific interactions and apply the formalism of the previous sections to calculate 
the BSF cross-sections, and the rates for de-excitation or decay of bound states, where relevant. 
We consider the interaction of two scalar particles (i) via a light scalar boson (Sec. 5.1), and (ii) 
via an Abelian gauge vector boson (Sec. 5.2). 

In the instantaneous approximation, these interactions are described in general by a Yukawa 
potential.(Of course, an unbroken gauge symmetry gives rise to a Coulomb potential.) A 
Yukawa potential admits bound state solutions if < a/i, where a is the fine-structure constant 
of the interaction. On the other hand, the radiative formation of bound states via emission of 
a force mediator is kinematically possible if < (a^ -|- for Vre\ < a - which is when 

the Sommerfeld effect renders bound-state formation efficient - this is a much stronger condition. 
Provided that this condition holds, the distortion of the wavefunctions due to the non-zero mediator 
mass, from their Coulomb limit, is expected to be negligible. For simplicity, we shall thus perform 
our computations in the Coulomb limit. 

As is well known, in the presence of an attractive Coulomb potential 

V{r) = --, a>0, (5.1) 

r 


there is a discrete spectrum and a continuous spectrum of energy eigenstates. The continuous 
spectrum corresponds to the two-particle states, and is characterised by a continuous quantum 
number that stands for the expectation value of the momentum of the reduced system, k = /rVrei, 
with Vrei being the expectation value of the relative velocity. The discrete spectrum corresponds 
to the bound states, and is characterised by the integer-valued quantum numbers {n£m}. In the 
discrete spectrum, the expectation value of the momentum of the reduced system is At/n, with 
K = fia being the Bohr momentum. As we shall see, the parameter that essentially determines the 
efficiency of BSF is the ratio of the momentum expectation values of the bound and the scattering 
states. 



a 

^rel 


(5.2) 


The energies of the states of the discrete and the continuous spectra are 


Wp.n — rn --h En , 

2m 


WK,k = rn 


K2 ^ 


(5.3) 


^^For a classification of the low-energy effective potentials generated by long-range interactions, and a systematic 
renormalisation procedure of singular potentials, see Ref. [75]. 
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where P, K are the momenta of the CM of the bound and the two-particle states respectively, and 




The wavefunctions are given in appendix F.l. 


k2 ^ A«4 i 

^ 2^i 2 


(5.4) 


Useful integrals 

For the calculation of the amplitudes Atk-s-n and we will find it useful to define the integrals 

5i(q,p; Ar,P) = J ^S{q-K) J ‘^Si{r]iP+p) {27r)S{q°-p°-rj 2 P°) , (5.5) 

E: 2 (q,p; iF,P) = J ^S{q-K) J ^^ 2 ( 772 ^ - p) (27r)d(9° - / -t 77iP°) . (5.6) 

We evaluate Si and S 2 in appendix E. We will also need the following integrals involving the initial 
and final state wavefunctions 


21k,n(b) = J 

f (^27r)3 4(P + b) , 

(5.7) 

J’k,n(b) = J 

f P^"(p) '^k(p + b) , 

(5.8) 

Ak,n(b) = J 

f (^27r)3 P^ ‘^k(p + b), 

(5.9) 

3 

III 

f ^2,^)3 ^n'(P + b) , 

(5.10) 

s 

III 

f (-2^)3 P<(p)^n'(p + b) , 

(5.11) 

3^ 

3 

III 

f (-2^)3 P^^n(p) 7^n'(p + b) . 

(5.12) 


We evaluate I, JT” and 1C in appendix F, for the initial and final states of interest. We shall use the 
integrals (5.5) - (5.12) in sections 5.1 and 5.2. 

5.1 Scalar mediator 

We consider the interaction Lagrangians 

SCs,r = d^Xi + ^d^X2 d^X2 + - \rnlxl - ^^2X2 “ 

- - ^52^27^X2 > (5-13) 

5 Cs,c = ^^lX\ d'^Xi + df,x\ d^X 2 + - m?|xiP - TO2IX2P - 

-gimit^IXiP-52W2¥7|X2p • (5.14) 

In Eq. (5.13), xi ^-nd X 2 are real scalar fields, while in Eq. (5.14) they are complex, (p is a real 
scalar boson, with mass mi, m 2 , and 51,52 are dimensionless couplings. 

r^In Eqs. (5.13), (5.14) and (5.48), we omit the quartic couplings in the scalar potential, since they do not enter our 
calculations. 
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viQ +p ^ viQ +p' 


viQ+p viQ + p' 


W{p,p'-,Q) = 



V2Q-P mQ-p' 


mQ-p mQ-p' 


Figure 6. In the non-relativistic regime, the one-boson exchange is the dominant contribution to the 
Xi ~ X2 interaction. Consequently, the 4 -point function is approximated by an infinite ladder of one-boson 
exchange diagrams. The exchanged boson may be either scalar or vector. 


To lowest order, the interaction between xi X 2 is mediated by one-boson exchange, as 
shown in Fig. 6. Then 

igig2mim2 


W{p,p';Q) = - 


{p - p'Y - ml ' 


In the instantaneous approximation 


14 


w(|p-p'l) = 


»gig2Wim2 

(p-p')2-bm2 


From Eq. (2.85), we find the non-relativistic potential 

1 


V{v) = -r 


That is. 


iAmp J {2 tt)^ 

V{r) = - 




(fp 


iAmp, 


(27r)3 p2 -I- ml 


with a = 


5152 


(5.15) 


(5.16) 


(5.17) 


(5.18) 


r IGtt 

The interaction is attractive if 5152 > 0, i.e. it is always attractive between particles of the same 
species or particles and antiparticles, but it can be either attractive or repulsive between particles of 
different species. As already mentioned, for our computations, we shall consider the limit —>■ 0. 


5.1.1 Bound-state formation amplitudes 

The lowest order contribution to the perturbative part of the radiative BSF amplitude arises from 
the diagrams shown in Fig. 7. In our approximation, the entire BSF amplitude corresponds to the 
ladder diagrams of Fig. 8. Recalling Eqs. (3.16) and (3.17), the diagrams of Fig. 7 evaluate to 


{27:^5'^ {qi + q 2 - Pi - P 2 - P<^) 5l, 52 ) = 

- igimiSi{qi)Si{pi) { 2 tt)^S^{P^ + Pi - qi) ^2(52) { 2 tt)^S'^{p 2 - 52) 

- *52W2<§2(P2)<S2(52) + P2 - 52 ) Si{qi) {2Tr)‘^S‘^{pi - qi) . 


Equivalently 

+ P^V 2 P - P,mK + q,p 2 K -q) = -S{q;K)x 
[gimi Si{r]iP + p) {2TT fS‘^{q - p - r]2P^) + 52^2 S2{ri2P - P) (27r)^d‘‘(5 - p-b giP,^)] . 

Note that the leading-order correction to the approximation of Eq. (5.16), and similarly of Eq. (5.50), are of order 
(p® — . Due to the linearity of the Bethe-Salpeter Eqs. (2.53), (2.54), the corresponding corrections to the 

approximations of the Bethe-Salpeter wavefunctions, ^Q^ri{p) and $Q,q(p), that are introduced via Eqs. (2.68) and 
(2.69), would also be of the same order. As discussed in earlier sections, in the regime of interest, p^ ~ 0{p^/p), 
and such corrections would then be 0(p^/p^), i.e. of higher order than the corrections considered in Eq. (3.32) 
and which appear in the S-matrix elements (3.29) and (3.30). This confirms the consistency of the approximation. 
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Figure 7. The lowest order contribution to level transition amplitudes, including bound-state formation. 


Using the above, Eq. (3.31) can be expressed in terms of the ^ 1,^2 integrals defined in Eqs. (5.5) 
and (5.6), as follows 


•^trans(Q7P) — 

_ [gimi Si(q, p; K, P) (27r)^^3(q _ p _ + g 2 W 2 S 2 (q, p; K, P) (27r)3(5^(q - p -b 7?iPy)] 

5o(p;P)5o(q; K) 

~-2mp 1+ [gi(27r)^5^(q-p-g2Pv>)+g2(27r)^5^(q-p-b77iPc^)] , (5.19) 


where in the second step, we used the non-relativistic approximations of ,:ii, .::;27 given in Eqs. (E.23), 
(E.24). Inserting this into Eq. (3.29), we obtain 


M 


k— 




(Pp 

(^ 


1-b 


2mp 


'0n(p) gi<^k(p + g2Pv>) + g2 0k(p-giPv:.) 


(5.20) 


The p^ term in the square brackets, will prove to be important in the cases of identical particles and 
particle-antiparticle pairs (see below). In terms of the integrals (5.7) - (5.9), Eq. (5.20) becomes 




qiP]^^n{j]2P^p') ~t“ g2Tk,n( giPc/?) “b 


gi^k,n(g2P(/7) T g2^k,n( ^iPip) 

2mp 


. (5.21) 


Capture in the ground state 

For capture in the {100} state, Eq. (5.21) becomes 

[gi2ik,{ioo}(g2Pvj) + g2Tk,{ioo}(—giP^)] 

^ gi^k(g2Py)+ g2^k(-giPy) 

2mp 

where we used Eq. (F.22). We shall drop the correction in the coefficient of the Ik,{ 100 } integrals. 
The mediator momentum is |P,^| = fk — '^^{loo} = (1 + C^)^^/(2g) (c.f. Eqs. (3.40), (5.4)). Using 
Eqs. (F.39) and (F.41) for the integrals Ik,{ 100 } and ^k, we find 


■^k-)-{ioo} — —m^/^ I 


b^k-j-tioo} — —P-iC) 


12m? 
pk 


X [( gig2 - 92Pi)cos0 


+ ^ {(ffi^2 + 92 V 1 ) [(—1 + *C) + 2(2 — i() cos^ 9 ] + ^(51 + 52)(1 + *C)| 


(5.22) 
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Figure 8. The ladder diagrams giving the dominant contribution to bound-state formation and other level 
transitions. The mediator can be either a scalar or a vector boson. 


where 0 is the angle between k and P,^, and TZ{C,) is given in Eq. (F.35). (We emphasise that 
the above expression is not a consistent expansion in a, but rather only in We discern the 

following cases: 


• For a particle-antiparticle pair, or for identical particles, ffi = 52 = g, ? 7 i = t ?2 = 1/2 and 
g, = to/ 4; the first term in Eq. (5.22) vanishes, and we obtain 

A4k^{100} - 7^(C) Urel [< + (2 “ <) COS^ &] 

^ [(1-f O Po(cos 6 >)-h (2 - O P 2 (cos 6 >)] , (5.23) 

where in the second line, we decomposed the amplitude in partial waves. 

• For non-degenerate particles, the first term in Eq. (5.22) dominates, and 


/to\ 

' {gim - g 2 mf 

VM/ 

10710 


-, 1/2 


Pi (cos 9). 


(The factor inside the square brackets becomes equal to 1 in the limit gi = 92 , Vi ^ V 2 


(5.24) 

•) 


Capture in excited state with non-zero angular momentum 

Since for a pair of degenerate particles, the cross-section for radiative capture to the ground state 
is either or suppressed (as seen by comparing Eqs. (5.23) and (5.24)), we shall now calculate 
the amplitude for capture in the {210} state. In this case, — £{ 210 } = (4 + C^)k'^/{8g). 

From Eq. (5.21), and keeping only the leading order terms of Eq. (F.42) for the Pk,{ 2 io} integrals, 
we find 


■^k^.{210} — —ms/^ [ffl21kT210}(^2P<p) + '?221k.{210}(~''7lP¥>)] 


/ 51»?2 +52?7i/ 


/*C + 2\ 

V V1071O / 

' 3 (4 +(2)3 ' 

UC- 2 y' 


X [(fC -I- 2)Po(cos0) -I- 8(*C - l)P2(cos0)]. 


The J = 0 and J = 2 contributions to the above yield the squared amplitudes 


|-^k^.{2io},J=oP 

|■^k^.{ 210 },J= 2 p 


TO 


{gig2 + g2gif 

10710 

{giV2+g2m)^ 

10710! 


2 I 9 .JJ .5 ^9 g—4Carccot((;/2) 

32 (4-f C^)4 1 - e-2^C ’ 

22^71^ ^9(1 -I-^2) g-4Carccot(C/2) 

32 52 (4-bC^)5 1 - e-2< 


(5.25) 

(5.20) 


(We remind that Alj is defined in Eq. (3.44), and note that the factor inside the square brackets 
becomes equal to 1 in the limit gi = 32 , independently of ??i,? 72 -) 
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5.1.2 Bound-state formation cross-sections and partial-wave nnitarity 

Combining Eq. (3.46) for the partial-wave cross-section and the amplitudes (5.23) - (5.26), we 
obtain 


• For gi= g 2 = g and = 772 = 1/2 (i.e. /i = m/4), 


G 


{ 100 } 
BSF, J = 0 


^{ 100 } 
BSF, J = 2 


crllOO} ^ { 100 } 

BSF BSF,J=0 ^ BSF,J = 2 


2^7r^ g-4CarccotC 

3^/i^ 1 -I- 1 ~ ’ 

2^7r2a3 ^4(4 ^2) g-4CarccotC 

325^12 (1 -^(2)2 1 _ e-2’^C ’ 

2^7r2a3 C'‘(3-H2C^) 

Ibfj? (1 -I- ('2)2 1 — e“2’rC 


(5.27) 

(5.28) 

(5.29) 


• For non-degenerate particles, 


grlioo} ~ 

BSF BSF,J = 1 


2 ^ 7 r 2 a 


{ 91 m - 92mY 


e e 


—4^arccot(^ 


• For capture to the {210} state, 
d 210 } ^ 


a 


BSF, J=0 


a 


{ 210 } 


gr{ 210 } ^ { 210 } 

BSF BSF,J = 0 


,.{ 210 } 


3 /r 2 

1677a 

(l-hC ^)2 l-e- 2 ’^C • 

28772 a 

'{9192 +9291)^' 

^8 g— 4 Carccot((^/ 2 ) 

32 /r 2 

214,72 a 

1677a 

’(51^2 +ff 2 ??l)^ 

(4 + C 2)3 1 — e“ 27 rC ’ 

^8^2^ _|_ (^ 2 ^ g- 4 Carccot(C/ 2 ) 

325^72 

28772 a 

1677a 

{gi 92 +9291)^ 

( 4 -KC 2)4 l-e- 2 ^C ’ 

C®( 28 -h 23 C^) e“ 4 Carccot(C/ 2 ) 

15/42 

1677a 

( 4 -hC ^)4 l-e- 2 ’"C 


(5.30) 

(5.31) 

(5.32) 
.(5.33) 


Unitarity limit 

For C > 1, Eqs. (5.27), (5.28), (5.30), (5.31) and (5.32) become 
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BSF, J = 0 
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BSF, J = 2 
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2'^Tra^ 
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(degenerate particles: gi = 92,91 =92 = 1 / 2 ) , 
(degenerate particles: gi = 92,91 =92 = 1 / 2 ) 
(non-degenerate particles) , 


{ 9192 - 9291 ) 


167ra 

(5i?? 2 +52»7i)^'' 


167ra 


2^‘^7ra^ \{9i92+g29if 


325268 


IGna 


(5.34) 

(5.35) 

(5.36) 

(5.37) 

(5.38) 


It is interesting to note that in this low-velocity regime, the velocity dependence of all partial waves 
is the same . ^8 This is in fact expected by unitarity, since 1 is both the large coupling and the 
low-velocity limit. Indeed, the unitarity bounds on the partial-wave inelastic cross-sections, shown 
in Eq. (3.47), all have the same velocity dependence. They are realised when the factors to the 
right of the x symbols in Eqs. (5.34) - (5.38) become « 1. The validity of our calculation is thus 
limited to at most a < Ouni) with the strongest bound, aun, « 1, obtained from CT,^gp°j_g- 

^®This was also noted in Ref. [66], for Sommerfeld-enhanced annihilation processes. 
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Figure 9. Annihilation of an unbound particle-antiparticle pair or decay of a bound particle-antiparticle 
pair into two force mediators. The mediator may be either a scalar or a vector boson. 


5.1.3 De-excitation rate 

The radiative capture to the {210} state is the dominant BSF process for particle-antiparticle 
pairs and pairs of self-conjugate particles. Moreover, for non-degenerate particles, it is slower but 
comparable to the capture to the ground state. Here, we shall thus compute the de-excitation rate 
of the {210} state. 

The radiative de-excitation of a bound state arises from the same diagrams as the radiative 
capture to a bound state, albeit for different initial and final states. In our approximation, these 
are the ladder diagrams shown in Fig. 8. Inserting Eq. (5.19) into Eq. (3.30), we find 


Mr 


~ -m { giTn',n{V 2 P ip) + g 2 ln',ni-'r]lPp) + 


gilCn',n{V 2 Pp) + g2JCn'.nj-m'P v) 
2mfi 


. (5.39) 


In the {210} —)■ {100} transition, the mediator is emitted with momentum |P,^| = £’{ 210 } ~ 
^lioo} = Then, using Eq. (F.45), we find 

-^{210}->'{100} — [5l2i{210},{100}(^2Pvp) + 1722i{210},{100}(~^lPv3)] 


i 2®m 

34 


gim + g 2 m 


•\/l67ra 

The de-excitation rate, given by Eq. (3.48), becomes 


r{210}-J-{100}-l-¥3 — 


2'^a^u r {gig2 + 52?7i)^’' 


37 


Wira 


(5.40) 


(5.41) 


5.1.4 Annihilation vs bound-state formation for particle-antiparticle pairs 

A particle-antiparticle pair or a pair of self-conjugate identical particles coupled to a light scalar, 
can annihilate into two mediators, x^X or XX shown in Fig. 9. In both cases, the 

perturbative annihilation amplitude, evaluated in the CM frame, in the non-relativistic regime, is 

■^ann(q:-q;P)-p) Ci ao -b -^(q- p) (5.42) 

mg 

with oo — oi ~ IOtto. Using Eqs. (4.24), (4.26) and (4.27), we find that the annihilation cross- 
section is s-wave dominated, o’ann'Crei — erg S'o.ann, with (To = /s|aoP/(327rm/j,) ~ jg? being the 

perturbative s-wave annihilation cross-section, and S'o,ann = |<('k(0)P being the s-wave Sommerfeld- 
enhancement factor for annihilation. Using the wavefunction of Eq. (F.ll), we obtain 

O’ann'Crel “ ^2 ^ ' (5.43) 


Note that here, g = m^l2 = m/4, with being the mass of the annihilating particles. 
We may now compare the BSF and annihilation cross-sections. 


cri^oo} 

BSF 

2®a^C^(3 -b 2/'^) 4farccotC 

(5.44) 

<^ann 

“ 15(1 + C2)2 

^{210} 

BSF 

^ 27C6(28 + 23C2) _4Carccot(C/2) 

(5.45) 

^ann 

15(4-bC^)4 
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At C > 1, 0-{i°°>/(Tann ^ End /(Ja.nn ^ 0.066d® We Compare tXann and cr|sp°^ in Fig. 1. 

Since annihilation is the dominant inelastic process, it lowers the value of a at which the unitarity 
bound appears to be realised, to auni ~ 0.54. 


5.1.5 Particle-antiparticle bound-state decay rates 

From (4.25), (4.26) and (4.28), we find that the decay of the {100} particle-antiparticle bound state 
into two mediators, is dominated by the s-wave contribution, 

'KO? 

r{ioo}^w = crolV'ioo(0)P = . (5.46) 

For the {210} bound state, the s-wave decay mode vanishes, since ^” 210 ( 0 ) = 0. However, a non¬ 
vanishing contribution arises from thep-wave mode. From Eq. (4.26), we find ui = /s|aip/(2^37rm/r) 
7ra^/(12p,^). From Eq. (4.28) and the wavefunction (F.9), we obtain S'l^dec = K®/(327r/r^). Then 

r - c _ TTO^ 

{210}^V>V l,dec 12/^2 327 rp 2 2^3’ ’ 

The decay rates into three mediators are expected to be suppressed by one additional power 
of a with respect to the above. Recalling Eq. (5.41), this suggests that for the excited state {210}, 
the transition to the ground state is the dominant decay mode. 


5.2 Vector mediator 


We now consider two scalar particles yi, X 2 coupled to a gauged U{1) force. The interaction 
Lagrangian is 

6Cv = {Dfj,xi)HD^Xi) + {Df,X 2 )HD^X 2 ) - , (5.48) 

where yi, X 2 are complex scalar bosons, F^'^ = — d^ip^ and — icjgip^, with ci, C 2 

being the charges of yi, y 2 -^^ 

The one-boson exchange diagram gives 

r.2 


W{p,p'-,Q) = 


*CiC25* 


[{2piQ +p + p')- (2772(5 -P- p')] ■ 


{p - p'Y 

In the non-relativistic regime, we shall approximate the above with 

urf I ciC2g‘^ 

W(p,p ;Q) ~-. 

(p - p')2 

In the instantaneous approximation, and setting ~ 

*4pi?72 CiC25^ 


W(lp-p'l):^-^ 


(P-P')^ 


(5.49) 


(5.50) 


and 


That is, 


V(r) 
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iAmp 


<fp 

(27r)3 


W(p)e*P'- 


= CiC25^ 


(fip 

(27r)3 p2 


V(r) 


a 

•> 

r 


with 


a = — 


ciC2g^ 

47r 


(5.51) 

(5.52) 


The interaction is attractive if C 1 C 2 < 0. 

^^Note that for a Dirac fermion-antifermion pair, the annihilation into scalars is dominantly p-wave. The spin- 
averaged annihilation cross-section times relative velocity is (Tann'^rei — cri‘S'i,anin with cri ~ 37rQ^/(2p^) and 
5'i,ann = C^)27r^/(1 — On the other hand, the BSF cross-sections do not depend on the spins 

of the interacting particles, since in the non-relativistic regime, the spin is conserved separately from the orbital 
angular momentum. Thus, for Dirac fermions, at ^ 1, 'TBSP^V'T'ann 0.1 and (T^s}°V'rann ^ 0.044/«2. This 
means that for a < 0.2, BSF is faster than annihilation in the regime where the Sommerfeld effect is important. 
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5.2.1 Bound-state formation amplitude, cross-section and partial-wave nnitarity 

The perturbative part of the level-transition amplitudes is with the lowest 

order contribution depicted in Fig. 7, 

{2Tr)"^6'^{qi + q2 - Pi - P2 - P^)iC^_^^p{P^,pi,P2,qi,q2) ^ 

- igciiPi + Qi) Si{pi)Si{qi) {2TTfd‘^{P^ +pi - qi) ^ 2 ( 92 ) {27r)'^6'^{p2 - 92 ) 

- igc2{p^ + ^ 2 ) S2{p2)S2{q2) +P 2 - ^ 2 ) Si{qi) (27r)‘‘5'‘(pi - qi) . 


From this, we obtain 


C^-amp(-P¥>) VlP + g 2 P “ P, VlK -t q, P 2 K - ?) = 

= -g S{q]K) |ci[2??iiF^ - (ryi - t] 2 )PI^ + 2p>^] Si[piP + p) {2TT)'^5^{q - p -'q 2 Pv) 

+ C2[2iq2K^ + {vi - g2)Pi^ - ‘^■P^]S2{'n2P - p) - P + ViPip)] ■ 

Then, Altrans = eM-^o-ans; where for j = 1, 2,3, 

■A^L„s(q;P) = - s,{ci;KUo{p-,P) 

X {ci [2771 - (771 - p2)Pi^ + 27F] E:i(q,p;iF,P) ( 27 r) 353 (q- p - 772?,^) 

-h C2 [2772-h (771 - 772)- 2 ^-^] S 2 (q,p;iF,P)( 27 r) 3 ^ 3 (q-p-H 77 iP<^)} . ( 5 . 53 ) 

We remind that Si,S 2 are defined in Eqs. (5.5), (5.6). Using their non-relativistic approximations 
(E.23), (E.24), we find the amplitude of Eq. (3.29) to be Atk->.n = where for j = 1, 2, 3, 

= -gi/^ J 7/;*(p) [2771 - (771 - r]2)P^ + 2 p>] 0k(p + V2P^) 

+ — [2r]2K^ + (771 - 772)P^ - 2p^] 0k(p - biP^)} • 

^2 J 

We may rewrite the above in terms of the X, J integrals, defined in Eqs. (5.7) and (5.8), as follows 


MU, 


= -25V^ ^ ^ - I JU-ViP^) 

PU Un{V2P^)+C2 (k^ + 


Cl - 


m 

m - V 2 
2?7i 


vi - m 

2772 


P^ Pk.n(-r7iP^) 


. (5.54) 


Because the vector boson Lp^ is transverse, the p = 0 component and the component parallel to 
do not contribute to the amplitude A4k-in = Dropping those components, we obtain 

MUn MUn- frame, K = 0, and for capture in the ground state {100}, using 

Eq. (F.40) for the Jk.fioo} integrals, and keeping only the leading term, we find 


= (5.55) 

where the sin0 factor arises from the projection of k on the plane vertical to P<p, and P(C) is 
defined in Eq. (F.35). (Recall that ciC 2 < 0 for an attractive potential.) Note that the partial wave 
decomposition of the 0-dependent factor of Eq. (5.55) is 

- - irr ^TT 

sin0 = vl — cos^ ^ ~ 4 Po(cos0) — P2(cos0) —^ P4(cos0) -|-... . (5.56) 
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The sum over the vector-boson polarisations is 

e e 


(5.57) 


In Eq. (5.57), the contribution of the /r = 0 and v = 0 components cancels the contribution from 
the component of A4 that is parallel to yielding the final result. Using Eqs. (5.55) and (F.36), 
we obtain 






{mci - mcif 


-C 1 C 2 


(5.58) 


Note that for ci = —C 2 = 1, which includes the case of a particle-antiparticle pair, the factor in the 
square brackets in the above expression becomes [(77201 — ? 7 iC 2 )^/(—C 1 C 2 )] = 1 . 


Using Eqs. (3.42) and (5.58), we find the unpolarised cross-section for radiative capture to the 
ground state to be 
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{ 100 } 

BSF 


X^rel 


( 772 C 1 - ?7iC2)^ 
-C 1 C 2 


2®7r^ C,^ g-4CarccotC 

3712 (TT^ 1 - e- 2 < ■ 


(5.59) 


This result agrees with Ref. [63] (see Eq. (75.6)), and is smaller than that of Ref. [64] by a factor 
of 2 . At C 3> 1, Eq. (5.59) becomes 


(T 


{ 100 } 

BSF 


^rel 


(772C1 - 77102)^ 
-C 1 C 2 


2Vg,2 

36 "^ 772 


(5.60) 


The contribution from the J = 0 partial wave is crli)!®] „ = ( 7 r/ 4)2 
unitarity limit on the s-wave inelastic cross-section, given in Eq. (3.47), is 
with 


0.69 


-C 1 C 2 


.( 772 C 1 - 77102)2 _ 


1/3 


(c.f. Eq. (5.56)). The 
realised for a = auni, 

(5.61) 


5.2.2 Annihilation vs bound-state formation for particle-antiparticle pairs 

In the non-relativistic regime, the annihilation of a particle-antiparticle pair into two vector media¬ 
tors, XX^ ipip, is s-wave dominated, with the perturbative unpolarised cross-section times relative 
velocity being erg = From Eqs. (4.24) and (4.27), we find 


<^ann77rel — ^0 |0k(O)| 


7ra2 27rC 
2/72 I _ g-27rC ’ 


where here p = m^/2 = to/ 4, with to^ being the mass of the annihilating particles. 
We may now compare the BSF and annihilation cross-sections: 


^{ 100 } g—4C arccotC 

Cann 3(1 -f C^)^ 


(5.62) 


(5.63) 


While at (/ < 1 BSF is negligible, at ^ ^ l-H the contribution of BSF to the inelastic scattering 
dominates over annihilation. At 1, cr ,^sF°Vo'ann 2:^ 1.56.^^ We Compare a^nn and in 

Fig. 1. 

^^For a fermionic particle-antiparticle pair coupled to a gauged U{1), ao ~ I(Afj?). Then, at 1, the 

importance of BSF relative to annihilation is even greater, with (Tggp*^^/(Tann — 3.1. 
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5.2.3 Particle-antiparticle bound-state decay rate 

From Eq. (4.25), we find the unpolarised decay rate of a particle-antiparticle bound state into two 
vector mediators to be 

l; /n\l2 fj,a^ 

r{ioo}^¥>¥> — CTo|'0ioo(O)| ~ 2^^ ~ ' (5.64) 


6 Discussion 

The formation of bound states affects the phenomenology of dark matter in a variety of ways. 
Computing the rates for bound-state formation and other related processes is essential in calculating 
the cosmology of DM and accurately estimating the expected DM signals and detection prospects. 

In the non-relativistic regime, the formation of bound states is enhanced by the Sommerfeld 
effect. The Sommerfeld effect has already been incorporated in computations of the DM annihilation 
rate, in the context of various theories, and has been shown to have important phenomenological 
implications. Besides enhancing the total DM annihilation rate, it may also modify - depending on 
the nature of the DM interactions - the relative strength of the various annihilation channels, thus 
changing the spectrum of the annihilation products [73]. Our results demonstrate that, for particle- 
antiparticle pairs or pairs of self-conjugate particles, the radiative formation of bound states can 
be faster than annihilation, in the entire regime where the Sommerfeld effect is important. This 
suggests that bound-state formation and decay may affect the annihilation signals of symmetric 
thermal-relic dark matter, as well as its relic abundance, well beyond the experimental uncertainty 
in the DM density [43]. Bound-state dynamics should then be incorporated in any relevant analyses. 

The importance of this point is underscored by present experimental results, which strongly 
constrain sub-TeV DM with electroweak interactions, and thus motivate investigations in the multi- 
TeV regime. As is well known, for symmetric (or self-conjugate) thermal-relic DM heavier than 
a few TeV, including WIMP DM, the Sommerfeld effect is important, both in the determination 
of the relic abundance and in the estimation of the expected indirect-detection signals. For the 
indirect detection of hidden-sector DM, the Sommerfeld effect - and therefore the formation of 
bound states - can be important even for lower DM masses. 

Asymmetric dark matter can couple even more strongly to light force mediators than symmetric 
DM; indeed, in the presence of a particle-antiparticle asymmetry, the very efficient annihilation 
which such a coupling would imply, cannot destroy the DM relic abundance. It follows that, for a 
much larger range of masses, asymmetric DM may efficiently form stable bound states in the early 
universe. This has important implications for its phenomenology. On one hand, the formation 
of bound states typically curtails the DM self-interactions and hastens the kinetic decoupling of 
DM from dark radiation in the early universe; consequently, it regulates the potential effect of the 
DM dynamics on the galactic structure. On the other hand, DM may participate in a variety of 
radiative processes inside haloes, such as excitations and de-excitations of bound states, or outright 
formation of bound states. In addition, the scattering of DM on nucleons may involve a variety 
of interactions, including both elastic and inelastic processes. This interplay between cosmology 
and the fundamental interactions of the DM constituents, determines all manifestations of DM 
today, and can be calculated only with precise knowledge of the rates governing bound-state-related 
processes. 

In this work, we established a field-theoretic framework for computing rates for processes in¬ 
volving bound states. This framework can be employed in future investigations of related effects, 
in a variety of theories. In particular, the computation of bound-state formation rates in theories 
which involve non-Abelian interactions - including the electroweak theory of the Standard Model 
- necessitates adopting a field-theoretic formalism. Moreover, this framework allows for systematic 
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expansions in the interaction strength and in the momentum exchange between the interacting de¬ 
grees of freedom; these higher-order corrections are important when the leading-order terms cancel, 
as was explicitly shown in our computations. 

The significance of long-range interactions - and therefore, the importance of comprehending 
their implications - is affirmed by unitarity. Unitarity sets an upper bound on the partial-wave 
inelastic cross-section, (crinei,j)max, shown in Eq. (3.47). This, in turn, yields an upper bound on the 
mass of thermal-relic DM [43, 71]. Notably, the velocity dependence of (crinei.j)max suggests that the 
unitarity bound can be realised only if the underlying interactions are long-ranged [43]. However, 
in the presence of long-range interactions, the formation of bound states may be the dominant 
inelastic process, as shown in the present work. The realisation of the unitarity bound, and its 
phenomenological implications, are thus largely determined by the dynamics of bound states, which 
should be fully incorporated in any related study. For example, the DM self-destruction in the early 
universe via bound-state formation and decay involves an interplay between capture, disassociation 
and decay processes that is absent in the case of direct annihilation into radiation [43]. 

Moreover, our results show that in the large-coupling (or low-velocity) regime, the dominant 
inelastic channel often belongs to a higher partial wave than usually assumed. Since the unitar¬ 
ity bound on higher partial waves is more relaxed, this implies that thermal-relic DM may be 
significantly heavier than previously estimated. 

Lastly, using the computed bound-state-formation cross-sections, we may estimate the interac¬ 
tion strength for which the unitarity bound is seemingly realised. Our leading-order computations 
show that this is at a ^ 0.5, i.e. well below what is often considered to be the perturbativity limit, 
a ^ TT or 47r. 
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Appendices 

A Bethe-Salpeter wavefunctions: Summary of definitions 


We summarise the definitions of the Bethe-Salpeter wavefunctions and their Fourier transforms. 
Following Sec. 2.1, we find 


^'Q,«(a7l,X2) 

= e *'^^4'q,„(37) , 

(A.l) 

^q.„ ( 3 ^ 1 . 272 ) 

= e^'^^^^q.Jo;) , 

(A.2) 

^Q,q(07l, 372) 

= e-*'5^$q,q(37) , 

(A.3) 

^q,q(a7l, 072) 

= e*'5^$*Q,q(:r) , 

(A.4) 


with = w„(Q) and = wq.q being the energies of and and where we defined 


^'q,«(a;) = (rj|rxi(?72a;)x2(-??ia:)|Sq,n) , (A.5) 

^q.n(a;) = (^Q.n|'7"xl(’72a;)x2(-’7ia;)|tl) , (A.6) 

^Q.q(a;) = {^\Txi{mx)X2{-'nix)\UQ.^) , (A.7) 

^Q.n(a;) = {UaJTx\{ri2x)xl{-mx)\^) ■ (A.8) 


The Fourier transforms are 


^Q.n(p) = ^ 

/’dSTq,„(3:)e*P" . 

(A.9) 

^Q,niP) = ^ 


(A.IO) 

®Q,q(p) = ^ 

[ d^37$q,q(a:)e*J>", 

(A.11) 

®Q.qW = j 

f dS<f>^,q(7r)e-*^>y 

(A.12) 


B Calculation of h{x^) 

As in Eq. (2.9), we define 

xi = X + r] 2 X , X 2 = X — r]ix , 

yi = Y + r]2y , y2 = Y -yiy . 

Then 


min(a;'j’,a;2) - max(j/°,?/2) = 

= min(X° -t r] 2 X°, X° - r]ix°) - max(r° -h 772 ?/°, Y° - yiy°) 

= X° + min(r72a7°, —riix°) — Y^ — max{r]2y^, —Viy°) 
f (-771 - 7721/°) , ifa;°>0,7/°>0 

(-77137° + 7?i7/°) , ifa;°>0,7/°<0 

(77207° - 7727/°) , ifa7°<0,7/°>0 

(77207° +7?iy°), ifa7°<0,7/°<0 

= X° - y° - i [(771 - 772 ) 37 ° -h (771 -f 772)|a7°|] -h i [(771 - 772)y° - (771 -h 772)|y°|] 
= X° -Y° + h{x°) - h{y°) , 


= x° - y° < 


(B.l) 

(B.2) 


(B.3) 
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where, as defined in Eq. (2.36) 


h±{x°) = i(?72 - Vi)x° ± . 


(B.4) 


C Calculation of iSo(p;P), iS(t;p;P) and iSi,2(^;p;P) 

Here we calculate the integrals 


5o(p;P) = j 

f ^ Sii-qiP+ p)S2ir]2P - p) , 

(C.l) 

S{t;p;P) = j 

f Si{r]iP+ p)S2{v2P - P)e"’'°\ 

(C.2) 

Si{t;p;P) = J 

f ^MmP + p)e-^^°\ 

(C.3) 

S2{t;p;P) = J 

f ^S2{v2P-p)e-^^°^. 

(C.4) 


(Note the different sign in the exponential between the definitions (C.2) and (C.3), (C.4), chosen 
so for later convenience.) 

We shall use the perturbative propagators 


Siipi) 


I 

2 2 i ^ ’ 

— mf + le 


S2(P2) = 


i 

~2 2 i ^ ’ 


(C.5) 


and we define Ei(p; P) = and E2(p; P) = -\/(772P^^li)M-Pr|. For convenience, 

we also define 


pi(p;-P) = - £’i(p;P), (C.6) 

P2 (p;P) = 772P°-£^2(p;P), (C.7) 

cri(p;P) = ?7iP°+ £'i(p;P), (C.8) 

cr2(p;P) = 772^°+ E2 (p;P) . (C.9) 


We first consider the integral (C.2), 


5(<;p;P)= f %- 


27r [?7iP0 +p° - Ei(p; P) + ie] [771 P° +p° + Ei(p; P) - ie] 

^ pip°t _ 

[772P° - + £'2 (p; P) - fe] [?72P° - - P2(p; P) + *e] 


(C.IO) 


We may evaluate 5(t;p;P) by closing the contour above and below the real axis for t > 0 and 
t < 0, respectively. The poles that contribute in each case are 


t>0: p° = -77iP°-Pi(p;P)+7e, = r72P° - P2(p; P) + 7e , (C.ll) 

t<0: p° = -77iP° + Pi(p;P)-ie, = r 72 P° + P 2 (p; P) - *£ • (C.12) 


Then 


5(t;p;P) 


g-ipi(p;P)t gicr2(p;P)t 

2 Pi [(po - Pi )2 - P|] + 2P2 [(p° + E2Y - El]\ ’ 

g-io-i(p;P)t giP2(p;P)t 

2 Pi [(po + Pi )2 - p2] + 2£;2 [(po _ E2Y - El]\ ’ 


for t < 0 , 


for t > 0 , 


(C.13) 
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and 


(C.14) 


So{P',P) = S{t = 0;p;P) = 


i(Ei + E 2 ) 


2EiE2 [(P0)2 - {E, + E2)^] 
It will be useful to rewrite Eq. (C.13), using Eq. (C.14), as follows 


/P° + Ei+E 2 
\P0-Ei+E2 
f P^-Ei-E2 

\P^ + Ei-E2 






P^-E^- E 2 
PO - El + E2 
P° + El + E 2 
PO + El - E2 


El 


El 


for t < 0, 


for t > 0. 


(C.15) 


(Note that in (C.13) - (C.15), it is implied that Ei and E 2 come with the arguments (p;P).) 
Similarly to the above, we find 


5i(t;p;P) = 


and 


52(1; p;P) = 


gio-i(p;P)t 

for 

t < 0 

2Ei(p;P) ’ 

gipi(p;P)t 

for 

t > 0 

2Ei(p;P) ’ 

g-ip2(p;P)t 

for 

t < 0 

2 E2 (p;P) ’ 

g—i(T2(p;P)t 

for 

t > 0 

2 E2 (p;P) ’ 


(C.16) 


(C.17) 


Non-relativistic approximation 

In the non-relativistic regime, P,p <C P°,mi,m 2 . Then 


P2\ P • p 


Ei(p;P) ~ ryi 
E2(p; P) ~ ry2 - 

Ei(p;P) + E2(p;P)-m + ^ + |^ 


m 

p p 

m 


2mi 

_p1 

21712 


(C.18) 

(C.19) 

(C.20) 


Note that in the last expression, the cancellation of the mixed terms proportional to p • P, can be 
traced to Eq. (2.11). This reflects the fact that in the non-relativistic regime, the relative motion 
can be separated from the motion of the CM. For convenience, we set 


pO = m+ — +£ . 


2m 


With these approximations, Eq. (C.14) becomes 

1 


5o(p;P) ~ - 


iimfi (P° - m - g) i4mn (s - g) 


(C.21) 


(C.22) 


D Partial-wave analysis for (co-)annihilation and decay processes 

We will prove the following two relations 


d^q 


(2^-M)!! 




dr^ 


dflr Pdcos6»p_r) <))k(r) 


(D.l) 


r=0 


- 39 - 
























and 


dQp Pii (cos Op) 


d?q 


( 27 r) 


3 </>k(q) |qrP£(cos6»q,p) = 

^ (2£+l)!! 


We shall use the addition theorem of spherical harmonics 


-P|(cos6»r) (^k(r) 


(D.2) 


J r=0 




m——£ 


(D.3) 


where x, y are unit vectors, and are the spherical harmonics. From Eq. (D.3) and the or¬ 
thonormality of Yirn, it follows that 

f 4:7T 

I (cos ^x,y) ~ 2 ^£ \ * (^•^) 

We will also need the expansion 


gjq r ^ -h l)i^je{qr) Pi^cosO^^r), 

i=0 


where jg is the spherical Bessel function, which satisfies 

d^ji{x) £l 

" (2^+1)!! ■ 


(D.5) 


(D.6) 


We begin with the right side of Eq. (D.l). Using Eqs. (D.3) - (D.6) and the Fourier transform 
of ^k(q)j we find 


d^ f 

^ £ / dDj. U^(cosdp^r) ((>k(r) 


47r 

2zn 

47r 

2ZTT 




J r=0 

d^q 


(27r 

d^q 


3 </>k(q) 


d^ 

dr^ 


dn^Yem{^r)P"^"‘ 


J r=0 


E f 7^<^k(q)E(2^' + i)*'' 

'i=-e ^ ’ e'=o 


/ 47r \ 


=-e 
2 i 




- r=0 


j E 


V2£-bl 
47ri^U f d?q 


d^q 

(27r)3 


<^k(q)|qn2^ + l)*' 


d^jejx) 

dx^ 


x—0 


0k(q) |q|^P£(cos6»p.q). 


(2£-bl)!!y (27r)3 
This proves Eq. (D.l). 

Acting on the left side of Eq. (D.l) with f dfip (cos dp), we find 

d?q 


I d£lj^Y£^(£lj^)P£f (cos dq r) 
Y^fn (^q) 

(D.7) 


dUp Pf'(cos dp) J 

(2£-hl)!! / 47r 


47ri^ £\ V 2£ + 1 

(2£-hl)!!/ 47r 
47ri^ £\ iv2£-b 1 ) 
, (2£-bl)!! 


{2£ -h 1) £! 


<('k(q) |qrPi;(cosdq,p) = 

J* dDp lY',o(^p) ^ £ J' dDj. P^(cos dp^r) 0k(p) 

i£ p 

/ dDr Y(!,o(f2r) (/'k(r) 
dDr P^(cosdr) <))k(p) 


J r=0 


r=0 


d*^ 

dr^ 


(D.8) 


J r=0 
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This proves Eq. (D.2). 


E Integrals for the non-relativistic reduction of transition amplitudes 


The transition amplitudes of Sec. 5 contain integrals of the following forms 

Ei{q,p;K,P) = J ^S{q;K) j ^Si{r]iP + p) {2Tr)6{q° - p° - r] 2 P°) , (E.l) 

S 2 (q, p;K,P) = S{q; K) J - p) (2^)5(g° - / + • (E.2) 


We are interested in evaluating Si(q, p-,K,P) at q — p — 'q 2 P^p = 0, and S 2 (q, p;K,P) at q — p + 
rjiPcp = 0. To evaluate Si and S 2 , we Fourier-transform the (5-function, and use Eqs. (C.15), (C.16) 
and (C.17). 

Eor Si, we obtain: 


Ei(q,p;iF,P)= dt J ^S{q-,K) J ^S,{r^,P+p) 

/ O /*oo 

(it iS(t; q; iF) iSi(t; p; P)-I- / dt S{t;q; K) Si{t;p; P) 

-00 J 0 


5o(q; K) 


where 

sr(q,p;E:,P) = 


[Pi(q;K) + P 2 (q;K)] 2Pi(p;P) 
i^o + Pi(q;K) + P2(q;K) 


X [Sr (q, p; K, P) + S+ (q, p; K, P)] , (E.3) 


P0-Pi(q;K) + P2(q;K) 
g°-Pi(q;K)-P 2 (q; K) 

iF0-Pi(q;K) + P2(q: 


P2(q;K) f dt 

J —00 


' —00 

/•o 


^Pi(q;K) f (E.4) 

7 J — 00 


and 


(q, p; K, P) 


iF0-Pi(q;K)-P2(q;K) 
iF0 + Pi(q;K)-P2(q;K) 
i^o + Pi(q;K) + P2(q;K) 
i^o + Pi(q;K)-P2(q;K) 


■P 2 (q;K) 

Pi(q;K) 



dt e®[“'^ 2 C”-l-p 2 (q;i<r)-l-pi(p;P)]t _ 


(E.5) 


Using the definitions (C.6) - (C.9) and the overall energy-momentum conservation, K = P + P^, 
we may simplify the phases appearing in the above integrals, as follows 

-772P"-pi(q;if) + (Ti(p;P) = P°-iF° + Pi(q;K)-fPi(p;P), (E.6) 

-d2^^° + 02(q;i^) + ai(p;P) = P" + P2(q; K) + Pi(p; P), (E.7) 

-772 P“-(7 i(q;P)+pi(p;P) = pO - 7^° - Pi(q; K) - Pi(p; P), (E.8) 

-V2P° + P2{ci;K)+p,{p;P) = pO-P2 (q;K)-Pi(p;P). (E.9) 


Adding a small imaginary parts to the integration variable t, such that the integrals in (E.4), (E.5) 
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converge, we find 


-i(q,p;if,P) i^ 0 _£;^(q.K) + £;2(q;K) 

/[i^o + i?i(q;K) + £; 2 (q;K)]ii; 2 (q;K) [K° - E^{<i;K) - E^iq^K)] E,{q;^) 

P° - KO + Ei{q;K) + E,{p;P) + E^iq^K) + E^ip-P) 

. [P° + E^ip-,P) + E^ (q; K) + i?2(q; K)] [E^ (q; K) + E^ (q; K)] - K^E^ (q; K) 
* [P° - KO + E,{q;^) + Ei{p;P)] [P° + E 2 {q-,K) + E^(p;P)] 

+ iq,p-, K,P) = j^o + E,{q-,K) - E2{q;K) 

- E^{q;K) - E2{q;K)] E2{q;K) [K° + E^;^) + E2{ci;K)] E,{q;K) 

P° - KO - Ei{q;K) - Ei{p;P) + P^ - E 2 {q;K) - E^ip-P) 

, [PO - i?i(p; P) - i?i(q; K) - £;2(q; K)] [i?i (q; K) + i?2(q; K)] - K°E, (q; K) 
■* [PO-ifO-Pi(q;K)-Pi(p;P)] [PO-P 2 (q;K)-Pi(p;P)] 


(E.IO) 


(E.ll) 


For So, we obtain: 


,(q,p;if,P)= ^S{q-,K) ^S2{mP-p) e-^P°*f{p) 

/ O roo 

dt S(t;q; K) S 2 {t;p; P) + / dt S{t;q; K) S 2 {t;p; P) 

-oo J 0 


5o(q;if) 


X [S 2 -(q,p;if,P) + S+(q,p;i^,P)] , (E.12) 


where 

(q, p; K, P) = 


[Pi(q;K) + P 2 (q;K)] 2^2 (p;P) 

+Pi(q;K) + P 2 (q;K) gdmP"-Pi(q;^)-P 2 (p;P)]t 

ifO-Pi(q;K) + P2(q;K) 

K — Pi(q, K) — P2(q, K) C ilrjiP°+a2{q;K)-p2{p-,P)]t CE IS") 

ifO-Pi(q;K) + P2(q;K)^'^'l’^V-oo ’ ^ ’ 


and 


-ti<l^P'^E,P) 


j^°-Pi(q;K)-P 2 (q;K) 
ifO + Pi(q;K)-P2(q;K) 
i^° + Pi(q;K) + P2(q;K) 
ifO + Pi(q;K)-P2(q;K) 


■i?2(q;K) 

Pi(q;K) 




dt e*[’'i-^^+^2(q;if)-CT2(pP)]t _ 


(E.14) 


Using the definitions (C.6) - (C.9) and the overall energy-momentum conservation, K = P + P^,, 
we simplify the phases appearing in the above integrals, as follows 


7;iP°-pi(q;i^)-p2(p;P) = -P° + Pi(q;K)+P 2 (p;P), (E.15) 

ryiP" + a2(q;/f)-p2(p;P) = -P° ++ P2(q; K) + P2(p; P), (E.16) 

77iP°-ai(q;i^)-a2(p;P) = -P° - Pi(q; K) - P 2 (p; P), (E.17) 

7;iP0+p2(q;i^)-a2(p;P) = -P° + - P2(q; K) - P2(p; P) ■ (E.18) 
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Adding a small imaginary parts to the integration variable t, as before, we find 


1 

^2 (q, p; P) = - ^0 _ (q. K) + (q; K) 

/[pO + Pi(q;K)+P2(q;K)]P2(q;K) - Pi(q; K) - ^^(q; K)] Pi(q; K) \ 

"" -P0 + Pi(q;K) + P 2 (p;P) -P0 + P0 + P 2 (q;K)+P 2 (p;P) ) 

^ . [P° - P 2 (p; P) - Pi(q; K) - P 2 (q; K)] [Pi(q; K) + P 2 (q; K)] - K) 

* [po-ifo i52(q;K)-P2(p;P)] [P°-Pi(q;K)-P2(p;P)] 

^+(q, p; K, P) = ]f^o + ^^(q.K)_£;2(q;-K) 

/[pO-Pi(q;K)-P 2 (q;K)]P 2 (q;K) + Pi(q; K) + P 2 (q; K)] Pi(q; K) \ 

"" ■ -P0-Pi(q;K)-P2'(p;P) -P0 + if0-P2(q;K)-P2(p;P) ) 

. [po + P 2 (p; P) + (q; K) + P 2 (q; K)] [Pi (q; K) + P 2 (q; K)] - /f0P2(q; K) 

* [p0-P0 + P 2 (q;K)+P 2 (p;P)] [P0 + Pi(q;K) + P 2 (p;P)] 

Non-relativistic approximation 


(E.19) 


(E.20) 


In the following, we consider the CM frame, K = 0, as in Sec. 5. We shall evaluate Si(q, p;P, P) 
and S 2 (q, p;P, P) at next-to-leading order in the momenta q, p, applying the non-relativistic ap¬ 
proximations of Eqs. (C.18) - (C.20), and setting, according to Eq. (C.21), 

= 771 + — -h , (E.21) 

2m 

= 771 + -h ^ (E.22) 


where and £\^ are given by Eqs. (F.7) and (F.13). The next-to-leading order corrections in 
q, p become important when the leading term in a expansion cancels, as is the case for the 
interaction of two degenerate scalar particles via a scalar mediator (c.f. Sec. 5.1). Note though 
that we drop subleading terms in the couplings; such corrections do not change the structure of the 
wavefunction convolution integrals (c.f. Eqs. (5.7), (5.8)) that enter into the transition amplitudes 
of Sec. 5, and thus do not avert the cancellation of the leading term in the expansion. The same 
holds for corrections. In addition, as seen from Eq. (3.40), P^ corrections are of order 
and while p^ corrections are only of order Similarly, P • p corrections are suppressed 

with respect to p^. As mentioned, we are interested, in particular, in evaluating Si(q, p;Ar, P) at 
q — p — r]2Pip = 0 and S2(q, p; AT, P) at q — p -I- "qiP^ = 0. Then, according to the above, we shall 
keep only the p^ corrections. 

At q — p — 772P,^ = 0, we find 


and at q 


-i(q,p; A:,P) ~ 2?7i2 5o(q; Ar)5o(p;P) 


p -I- r^iP^ = 0, we find 


■=‘2(q, p; K, P) ~ 27711 5o(q; K) 5o(p; P) 



(E.23) 


(E.24) 


F Schrodinger wavefunctions and convolution integrals 

In the following, we shall consider the attractive Coulomb potential 

F(r) = -“. (F.l) 

r 
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F.l Solutions of the Schrodinger equation 

The discrete spectrum of solutions to the Schrodinger equation (2.83), 




V’„(r) = Sn-ipnir) , 
for the Coulomb potential of Eq. (F.l), are (see e.g. [76]) 


Rneir) = 


(n — i— 1)! 
2n(n + i)l 


1/2 


{2nr/nY (2Kr/n) 


(F.2) 

(F.3) 

(F.4) 


where Yim{^) are the spherical harmonics and are the generalized Laguerre polynomials of 

degree n — i In the above, 

K = fjLa (F-6) 

is the Bohr momentum, and the energy eigenvalues are 

(F.7) 


f — — 


fia 


2v? ' 


In the following, we shall consider transitions to the ground state {100}. For the case of a scalar 
mediator, we shall also consider radiative capture to the first excited state with non-zero angular 
momentum {210}. The corresponding wavefuctions are 


■0100 (r) = 


V’2io(r) = 


k3/2 
!— ^ 

4v^ 


cos6lr , 


(F.8) 

(F.9) 


where 0r is the polar angle of the position vector r. 

The continuous spectrum of solutions to the Schrodinger equation (2.84) 

(()k(r) = £k</>k(r) , 




(F.IO) 


is characterized by the quantum number k, which is the expectation value of the momentum of the 
reduced system, with k = /rv^ei and Vrei being the expectation value of the relative velocity. The 
solutions are (see e.g. [76]) 


0k(r) = r(l — i^) F [/^, 1, i{kr — k • r)] 


where k = jkj, and 


( = n/k . 

2^JL~ 2 ■ 


The energy eigenvalues are 

Thus, C, = a/Urel- 

assume the following normalisation for the Laguerre polynomials 

r(n -h a -h 1) 


(F.ll) 

(F.12) 

(F.13) 


f x^e ^ Lm\x) dx = 

Jo 


(F.5) 
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F.2 Convolution of the wavefunctions 

We now want to calculate the integrals (5.7) - (5.9), appearing in the amplitudes of Sec. 5, 

(fp 

(27r)3 

cPp 


21k,n(b) = 
JT” k,n(b) = 

/Ck,n(b) = 


^n(p) <^k(p + b) = / d^r 


and 


21„',n(b) = 

Jn',n{h) = 
^n'.n(b) = 


(27r)' 

■ d^p 
(27r)3 

d?p 

(27r)3 

d^P 

(^ 

d^p 

(^ 


P V’n(p) <(>k(p + b) = i / ^3^ [VV'*(r)] (/)k(r)e 


-zb-] 


P V'n(p) '(>k(p + b) = - / d3r [V V’n(r)] <(>k(r)e 


<(p)^n'(P + b) = / d3r.i/;*(r)i/'„'(r)e 


ibr 


P i’niP) i’n'{p + h) =i / [V'0*(r)] ?/’„'(r)e 


—ibr 


P V'))(P) V’n'(p + h) = - d-^r [V V’n(r)] V'n'(r)e ' 


^ —ibr 


where we transformed into the coordinate space using Eqs. (2.81) and (2.82). 

For the {100} state, the integrals (F.14) - (F.16) become 
^3/2 girC/2 p^'J. — iQ 


2^k,{ioo}(b) — 
k,{100}(b) = — 


j^5/2 e</'^T{l-iC) 


j d^r e*(k-b).r-«r p - k • r)] , 

d^r f e*(k-b) r-«r ^ - k • r)] 


^k,{ioo}(b) — —K^21k,ra=i(b) +/Ck(b) , 
where f = r/r and 


^k(b) = 


2^5/2 g.C/2r(l_,^) 


F [j^, 1, i(kr — k • r)] . 


For the {210} state, we will need only the integral (F.14), 

^5/2 g7rC/2 r(l _ 


2^k,{210}(b) — 


4v^ 


To evaluate these expressions, we make use of the identity [64] 
r ^ ghk-b)-r-Kr [b2 + (^ - ifc)^] 

d r -^- FHA,^ikr-'k■r)] = ^ = /k,b(K) • 


Airr 


[(k-b)2 + K2]i-< 


(F.14) 

(F.15) 

(F.16) 

(F.17) 

(F.18) 

(F.19) 


(F.20) 

(F.21) 

(F.22) 

(F.23) 


J d^r r cos dr F [<, 1, i{kr - k • r)] . (F.24) 


(F.25) 


Differentiating Eq. (F.25) with respect to k [29], with respect to Id [64], and with respect to b, we 
obtain the following expressions 


Ik,{ 100 } (b) = r(l - iC) , 


'^k,{100} 




/Ck(b) = 8v^ r(l - iC) /k,b(K) . 


Ik,{ 210 } (b) = -tJ- r(l - zC) 


d'^fk,b{K') 
db Ok' 


(F.26) 

(F.27) 

(F.28) 

(F.29) 
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For the cases of interest, b = ri 2 P^ or —t]iP^-, evidently 6 = |b| < |P,p|. Moreover, |P,^| 
is determined by energy-momentum conservation (c.f. Eq. (3.40)). In the CM frame (K = 0), 
and in the non-relativistic regime, |P,^| ~ + ^rei)/^ for capture in the ground state, and 

|P,p| ~ /r(a^/4-|-u^gj)/2 for capture in the {210} state. As long as a, r^rei < 1, then & «C k if a > Urei, 
OT b k, if a < ?;rei- In evaluating Eqs. (F.26) - (F.29), we may thus expand in b and keep only the 
leading orders. In particular, to estimate 9/k,b(K)/9«;, dfk,h{K)/db^ and f\^_b[K’)/dbdK’, we first 
differentiate /, then use k = (k (c.f. Eq. (F.12)) and expand up to order 6^ around b — 0. (Note, 
though, that in most applications, we will not need all of the terms included in the expansions 
below.) Setting 

cos 6 = , (F.30) 

k b 

we find 


/k.b(K) 


1 

k^{i + C^y 



26 cos^ ^ 

1 1 1 

— (1 — iC) + 2(2 — iQ cos^ 9 

V<-i; 

k{i + iC) 

y{i + C^){i + iO 


(F.31) 


dh.h{K) 

4C6 (^<+1^ 


b 

^ 1 

— (1 — zC) + 2(2 — i() cos^ 9 

Ok 

k^{i++iC) v<-i/ 


k{i+e) 


(F.32) 


dfk,b{K) _ _2_ / iC + 1 

dW fc4(l-I-C^)(l + *C) \*C~1 


—K 


b> 

1 + fC 


1 + 


46 cos 9 

k{i+e) 


y 


2(2 - iC)bcosd 26^ Tl - - (6 - - ibQ cos^ 

k{i + c) + 


and 


’5^/k,b(K') 

2«K«-2) 

(^C + 2\ 

dbdK' 


{iC-2) 


— 2 — 4(i^ — 1) cos^ 9 — 
For convenience, we define 


86(f^ — 1) 3{i( — 2) — 4(z^ — 3) cos^ 9 


cos 9 


ki-i + C^ 


F(l-iC) + 


1-kC^ i-k< V*C-i 

In evaluating the cross-sections of interest, we shall need 

^6 p—4(^arccot 


|R(OI^= ' 


where we used the identities 


|F(l-iC)r = 27rC 


(l-kC^)3 ’ 

-ttC 


1 - e-2< ’ 


zC + 1 
zC-1 


-*C 


_ arccot 


(F.33) 
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Collecting the above, we obtain 


27^(C) h / ~ 6 -1 + zC + 2 (2 -Ocos2 6» 

Ik,{100} (b) [ cos 0 + —— 


(F.39) 
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46 cos 9 
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/Ck(b)~7^(C)fcl/2 


^ ^ 2b cos 9 

! + <+ — + 


62 


fc2(l + C2)2 


— 1 + iC + 2(2 — i() cos^ 0 

l?(TTc^ 


and 


i2'^V^b Cy^iC-2)e^</^T{l-tC) fiC + 2Y^‘^ 
k,{ 2 io}()^ p /2 (4 + C2)4 V<-2j 

8b{i( — 1) 3(iC — 2) — 4(i^ — 3) cos^ 9 
iC, —2 — 4:{i( — 1) cos^ 9 -^ 


(F.40) 


(F.41) 


cos 9 


ki4 + Y 


(F.42) 


Note that in the above, the angle 9 is related to the angle 9 between k and defined by 


as follows 


COS0 


k-P^ 

k\py ’ 


6 », for b = 772P,^ , 

TT + 0, for b = -ryiPcp . 


(F.43) 


(F.44) 


For the computation of the {210} state de-excitation rate (Sec. 5.1.3), we will need the integral 
1 (210},{ioo}(b). Starting from Eq. (F.17), and using the wavefunctions (F.8) and (F.9), it is easy 
to show that 


-r (26)/(3 k) i 2^72 6 

I{210},{100}(b) - cs 


(F.45) 


where we took into account that in the {210} —)• {100} radiative transition, the mediator is emitted 
with momentum |P,^| = {3/8)fia‘^; then, for 6 = 77i,2|P(p|, 26/(3«;) = t]i ^2 a/4 1. 
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